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ABSTRACT 


The goal of this investigation is to design minimum 
weight arch structures which span the distance between two 
points in two-dimensional space. An arch of unknown shape and 
variable cross-sectional width is modeled as a series of 
straight bar-beam elements. Finite Element Methods are used 
to compute the stresses in each element. Automated Design 
Synthesis (ADS) software is then used to vary the slope of 
each element and the cross-sectional width to prevent the 
yield stress of the material from being exceeded as ADS 
minimizes the arch volume to arrive at the minimum weight 
structure. Results are presented for a number of different 


loadings and boundary conditions. 
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I. INTRODUCTION 


A. MOTIVATION 

Since the beginning of recorded history, the arch has been 
used in the design of a wide range of structures. The arch 
was initially developed as a more reliable and efficient way 
of building a portal in a wall comprised of bricks or stone 
blocks. Over the years, it has been applied to an ever 
increasing number of load bearing structures. A significant 
reason for its popularity with architects and engineers is the 
aesthetically pleasing form of the arch. Many designers still 
enjoy the idea of an ordinary person being moved to wonder how 
a structure manages to stand with so little visible support. 

The development of the digital computer in the latter half 
of this century placed an incredibly powerful tool in the 
hands of designers who previously had to rely mainly on 
experience to select the best arch design for a given 
application. Using the high speed calculation capability of 
the computer, a designer could now use the techniques of 
numerical analysis to approximate the solution of a complex 
differential equation with a high degree of accuracy. These 
techniques could be used to model a structure under loading 
and obtain the optimum design. The arch has been the subject 


of several such studies. 


In 1976, Farshad [ref. 1] derived two governing nonlinear 
partial differential equations for statically determinate, 
hinged-hinged arches. The system’s total potential energy, 
modified by several objective functions through the use of 
Lagrange multipliers, was minimized in regards to design and 
state variables to achieve optimality and equilibrium. For an 
arch with specified span and loading, the nonlinear systems of 
equations for optimal thrust, minimum length, and minimum 
volume were presented, but not solved. 

In 1980, Rozvany et al. [{ref. 2] used the Prager-Shield 
criteria to determine the optimal shape of hinged-hinge 
frames. His frames consisted of two rigidly connected 
inclined beams with a point load applied at mid-span. He 
concluded that, for a single load condition, the optimal 
structure developed either bending only or axial forces only 
in the entire arch. When the ratio 4L/D, where L is the 
length of the span and D is the constant depth of the cross- 
section, is greater than eight, only axial forces develop. In 
this case, the optimal height of the arch is L/2. When 4L/D 
is less than eight, only bending develops and the optimal 
structure is a straight beam. The width of the cross-sections 
varied linearly from support to axis of symmetry. The author 
also found that although for a single system of point loads, 
the optimal arch consisted of straight segments, curved 
segments occurred if several alternative systems of point 


loads were considered. 


Also in 1980, Lipson et al. [ref. 3] conducted a numerical 
study using an automated design routine to determine an 
optimal arch design with the arch shape and cross-sectional 
dimensions allowed to vary. He modelled a uniformly loaded 
parabolic arch as a system of straight segments with thin- 
walled, rectangular cross-sections. When he studied arches of 
constant depth and width, the wall thickness determined the 
optimal shape for minimum weight. The resulting shape was 
found to be a parabolic curve with a height of 0.342 times the 
span length. 

In 1988, Ang et al. [ref.4] investigated the optimal shape 
of plastically designed non-funicular arches under a uniformly 
distributed load. The problem was solved by parameterizing 
the unspecified arch axis using spline functions and employing 
a smoothing function to approximate the non-smooth objective 
function (arch weight). The arch considered had a rectangular 
cross-section of variable width and constant depth. lis 
carried a uniform load and its supports were either hinged- 
hinged, hinged-clamped, or clamped-clamped. The optimum shape 
of the arch was found to be a parabola with a height of 0.433 
times the span length, which appears to disagree with Lipson’s 
results. 

In 1990, Charles Scott McDavid of the Naval Postgraduate 
School [ref. 5] investigated the optimization of circular 
arches subjected to various loading and boundary conditions. 


He modelled the arches as systems of straight segments with 


constant depth and variable width. He concluded that the 
bar/beam model was a viable technique in the approximation of 
arch structures and that the more statically indeterminate 

an arch structure was, the more efficient it was under 
identical loading. Lieutenant McDavid also proposed several 
further topics of research in this areas. One of the 
recommended topics was taken up Margaret Anne Menzies in 1991 
{ref. 6]. She investigated circular arches with varying depth 
and width. She also validated the bar/beam model as an 
approximation of an arch, but was limited by computer 
restrictions to a relatively small number of straight 
segments. This investigation looks into another area 
suggested by Lieutenant McDavid by allowing the radius of 
curvature (shape of centroidal axis) to be one of the design 


variables. 


B. PROBLEM STATEMENT 

The arch considered by this investigation has a rectangular 
cross-section with constant depth and varying thickness. The 
radius of curvature is allowed to vary to obtain the optimum 
shape of the arch. The distance spanned in the horizontal and 
vertical directions will be specified. The cross-sectional 
dimensions are assumed to be small with respect to the radius 
of curvature, which implies that the centroidal axis and the 


neutral axis coincide. 


The governing equations for the behavior of the arch are 
the beam equilibrium equation: 
(EIv")" = P(g) (1) 
and the bar equilibrium equation: 
(Agu!) = -P.(s) (2) 
The prime superscript indicates a derivative with respect to 


the independent variable,s. The variables in the equations 


are defined as: 


E = Young’s Modulus 

I = Cross-sectional Moment of Inertia 
A = Cross-sectional Area 

v = Lateral displacement 

u = Axial displacement 

P, = Lateral loading 

P. = Axial loading 


s = independent variable 


The arch is approximated by a system of straight elements. 
The local displacements are used to generate the internal 
pseudo stresses by applying the virtual load techniques 
described by Ding and Esping [ref. 7]. When the stress 
distribution is known, the volume of the arch is minimized to 
obtain a shape that keeps the developed stresses below the 


maximum allowable stress. 


The purpose of this study is to minimize the total weight 
of an arch spanning a specified horizontal and vertical 
distance under a variety of loading conditions. The material 
of the arch is assumed to be linearly elastic, homogeneous, 
and isotropic. The slope of the arch and the cross-sectional 
width at each node are allowed to vary to obtain the optimum 
(minimum weight) structure. Automated Design Synthesis 


software [ref. 8] is used to carry out the optimization. 


II. PROBLEM FORMULATION 


A. PROBLEM STATEMENT AND ASSUMPTIONS 

As stated previously, the goal of this investigation is to 
design a minimum weight arch which spans a_ pre-defined 
distance in two-dimensional space. The radius of curvature 
and cross-sectional width of the arch are allowed to vary to 
obtain the optimum shape. The following assumptions are made 
to simplify the problem: 

e The arch material is homogeneous, isotropic, and linearly 
elastic. This implies that the minimum volume arch is 
also the minimum weight arch. 


e The arch is modelled as a series of straight bar-beam 
elements whose behavior is governed by equations 1 and 2. 


@e Failure is assumed to occur when the yield strength of the 
material is exceeded. 


B. MATHEMATICAL MODEL 

The arch spans a distance in two-dimensional space from 
point A to point B (see fig. 2.1). The number of elements, 
NEL, the cross-sectional height, h, and the horizontal 
distance spanned by each element, DX, are specified. DX is 


constant and is defined as: 


DX = AX = L/NEL (3) 


where L is the distance between A and B in the horizontal 


direction. 





Figure 2.1: Arch Model 


The design variables in this case are the cross-sectional 
base, b,, at each node and the slopes, SLP,;, of the first NEL 


- 1 elements. SLP,; is defined as: 


SLP, = @, = (Y,,, - Y,)/DX 2 


where Y,,, and Y, are the vertical coordinates of the nodal 
points at the ends of the it® element. Since the structure 
must rise a distance H, the slope of the last element is not 
a design variable, but is fixed by whatever values the slopes 
of the previous NEL-1 elements take. 

The base dimension from one element to an adjacent element 


maintains smooth piece-wise continuity, as does a plot of the 


spatial coordinates of the nodal points. The resulting 
elemental shape is a three-dimensional trapezoid (fig. 2.2) 
whose volume is determined by multiplying the average base and 
the height with the length of the element. MThe elemental 


length, 1,, is defined as: 


1, = DX? + (SLP, * Dx)? (5) 


Therefore, the volume of the it® element is defined as: 


Volume; = D,,., * 2+ 1; 

where: 

Dave, = (Din. + D;)/2 

b, = nodal base dimension (6) 
h = nodal height dimension 

1, = elemental length 


i=i®™ element 





Figure 2.2: Elemental Shape 


The volume of the entire arch is obtained by summing the 


individual elemental volumes. 


C. OPTIMIZATION PROBLEM 

The goal of this investigation is to obtain a minimum 
weight (volume) arch which maintains a stress distribution 
which does not exceed the yield strength of the material. The 


volume of the arch is defined as: 
NEL NEL 
Volume = » Volume; = )° b,,., * H* 1; (7) 
sl fal 


This equation can be rewritten as: 


NEL 
Volume = h * DX * ()) dye, * 1 + SLP;?) (8) 
i=l 


The objective function for the optimization problem is defined 


ass 


NEL 
OBS = » Dayo, * YL (SEPA (9) 
sl 


Since h and DX are constant, minimizing the objective function 
will minimize the arch volume. 

Constraints are imposed upon the objective function and 
design variables to ensure the applicability of the governing 
equations and to prevent failure by yielding. The strength 


constraint can be simply stated as follows: 


10 


0,558, 


or in normalized form: 


y (10) 
where: 


0, = maximum stress at i“ node 


WY 
ul 


material yield strength 


The stress distribution in the arch is obtained by Finite 
Element Methods. This process is described in Chapter IV. 
The constraints placed on the base dimensions of the arch 
are used to ensure the validity of the bar and beam 
equilibrium equations (equations (1) and (2)). Limiting the 
base relative to the constant height prevents the structure 
from behaving like a shell ar a deep curved beam. The 
following limits are imposed through the use of the ADS side 


constraints on b; : 
Q.i*b sb, < 3*h (11) 


It is also important to choose a cross- sectional height 
dimension that is small relative to the length of the entire 
structure. 

One limit on the cross-sectional dimensions that was not 
considered by this investigation is that imposed by the 
requirement for elastic stability. The obvious importance of 
avoiding failure by buckling is an interesting area for future 


research. 


11 


III. OPTIMIZATION ANALYSIS 

The computer optimization is performed by ADS, a general 
purpose numerical optimization program containing a wide 
variety of algorithms. The program minimizes one function, 
the objective, subject to bounds imposed on the design 
variables by constraint functions. The solution of the 
general problem is separated into three basic levels which are 
strategy, optimizer, and one-dimensional search. For this 
constrained minimization problem, the ADS user’s guide [ref. 
8] recommends a Sequential Linear Programming strategy, a 
Modified Method of Feasible Directions optimizer, and a Golden 


Section Method one-dimensional search. 


A. STRATEGY 

Sequential Linear Programming (SLP), linearizes nonlinear 
objective and constraint functions, then obtains a solution to 
the approximation of the problem by linear programming methods 
(see fig.3.1). The problem is linearized again about the 
design point yielded by the solution to the previous 
approximation and resolved. This process is repeated until a 
precise solution is obtained [ref. 9]. 

The nonlinear functions are linearized through the use of 


a first-order Taylor Series expansion in the following manner: 


12 


Minimize: F(X) ~ F(X,) + VF(X,) * 6X 
Subject to: g,(X,) + Vg,(X,) * 8X < 0 as 
where: 6X=X- X, 

j = j constraint 
The zero subscript identifies the point about which the Taylor 
series expansion is performed. At the initial design, the 
objective and constraints are linearized to give straight line 

representations of the functions. 

In an under-constrained problem like this one, where there 
are fewer active constraints than design variables, this 
method sometimes performs poorly. This occurs because the 
linear approximation may be unbounded. This problem is dealt 
with by imposing move limits on the linear approximation as 
shown in figure 3.1. This ensures that the optimum is 
eventually reached within a tolerance of the move limits. In 
practice, these move limits are reduced during the design 
process so that a solution is found with the desired accuracy. 
SLP tends to converge rapidly to a solution, but while the 
solution of the linear problem is near the nonlinear optimum, 


it is in the infeasible region. 


B. OPTIMIZER 
The Modified Method of Feasible Directions for constrained 


minimization is used to solve the approximate optimization 


FS 


F = CONST, 


7 


INFEASIBLE REGION 
TRUE OPTIMUM 


APPROXIMATE OPTIMUM 


OBJECTIVE FUNCTION CONTOURS (F) 


var 6 CONSTRAINT FUNCTION 


LINEAR APPROXIMATION 





Figure 3.1: Sequential Linear Programming 


sub-problem (eq. 12). In this method, a search direction 
vector is first found. The vector containing the design 
variables is then updated by moving in the direction of the 


search vector in the following manner (refer to fig. 3.2): 


14 


NE ey, Oe, (13) 

The scalar a* defines the distance moved in the direction of 
the search vector in design variable space and q represents 
the iteration number. From the initial design, the objective 
function is moved in the direction of steepest descent in 
design variable space. This will continue until a constraint 
is encountered. Once this happens, a new search direction is 


obtained by solving the sub-problem: 


Maximize: -VF(X) * 3 
Subject to: Vg,(xX) * <0 jeg 
S+*8<1 
where: (14) 
X = vector whose elements are the design variables 


F = the objective function 

5 = the search vector 

J = number of active constraints 
The search direction will follow the constraint, but will 
allow the design to leave the constraint boundary if the 
objective function can be further reduced. If the scalar 
product of the gradient of each critical constraint with the 
search vector is less than zero, the search vector is moving 
away from the currently active constraint. This constraint is 
then dropped from the set of active constraints. If the 
search vector is the null vector, or numerically small, the 
optimization is terminated because this indicates that the 


Kuhn-Tucker conditions for optimality have been met [ref. 9]. 


1) 


| Hy TRUE OPTINUM 





Figure 3.2: Modified Method of Feasible Directions 


C. ONE~DIMENSIONAL SEARCH 

The scalar a* in equation (13) is found by the Golden 
Section Method. Basically, this determines how far the 
optimizer will search in the direction of the search vector. 


The one-dimensional search attempts to find the optimal value 
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of a* which will result in the minimum value of the objective 
function. It accomplishes this by progressively bracketing 
the minimum of the function by comparing the values of the 


function at smaller and smaller limits. 


D. ADS PROGRAM PARAMETERS 

Reference 8 contains a table of ADS program parameters. 
It is possible to change these parameters from their default 
values with the first call to ADS using an INFO value equal to 
-2. The purpose of modifying these parameters is to try to 
obtain a better optimal design by "fine tuning" the program. 
For instance, modifying the constraint tolerance, CT, can 
ensure that the active constraints will be closer to the 
desired value. However, this can also mean that it will take 
much longer for the solution to converge. The only parameter 
modification done in this investigation was to inhibit the 


auto-scaling function of ADS. 


di? 


IV. STRESS ANALYSIS 

In order to accomplish the goal of designing a minimum 
weight arch, a FORTRAN computer was used to evaluate the 
objective function and constraints and to invoke the ADS 
software. The strength constraints require that the stress at 
each nodal point remain below the yield strength of the 
Material. Determining the stress distribution in the arch is 
a complex problem, because of the need to be able to analyze 
a statically indeterminate structure. Finite Element Methods 
are used to calculate the displacements at each nodal point. 
The displacements are then related to the nodal forces and 
Moments by the direct stiffness method. Once the forces and 
moments are known, the stresses can be evaluated with a 


knowledge of the cross-sectional dimensions. 


A. STRESS DEVELOPMENT 

For the purposes of this investigation, the aren is 
assumed to fail if the total stress at any nodal point exceeds 
the yield stress of the material. The total stress is defined 
as the sum of the normal stresses due to bending and axial 


loading. 
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G, = |O,| + 19, 


where: 

0, = total normal stress (15) 
0, = normal stress due to bending 

0, = normal stress due to axial load 


The absolute values of the bending and axial stresses are used 
because at either the top or bottom extreme fiber the stresses 
will be additive. In other words, the yield stress is assumed 
to be the same in compression and tension. There will be 
shear stresses developed in the arch, but the geometric limits 
imposed by the side constraints ensure that these stresses 
remain negligible (see Appendix A). 

To calculate the two normal stress components, the arch is 
modelled as a series of straight elements, which can be 
considered to behave as beams to calculate the bending 
stresses and as bars to calculate the axial stress. The 


bending stress in an element is defined as: 


where: 

; 16 
M = bending moment =) 
c¢ = distance from centroidal axis to extreme fiber 


I = Moment of Inertia of cross-section 


The moment, M, is calculated by: 
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M = EIv" 
where: 
E 


Vv 


(17) 
Young's Modulus 


lateral displacement 

By substituting equation (17) into equation (16) the following 
equation is obtained: 

0, = Ecv" (18) 


In a similar manner, the axial stress is calculated. The 


stress due to axial loading is defined as: 


where: (19) 
F = axial load 
A = cross-sectional area 


The axial force, F, is determined by the following equation: 


F = AEu! 
where: (20) 
u = axial displacement 


By substituting equation (20) into equation (19), the 


following result is achieved: 


0, = Eu’ (21) 


Combining equations (18) and equation (21) results in: 


0, = E(cv" + u’) (22) 
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The approximate values of the lateral and axial 
displacements are obtained through the use of the Galerkin 
Finite Element Method. With the use of these values, the 


stresses at each nodal point can be calculated. 


B. FINITE ELEMENT BEAM EQUATION DEVELOPMENT 

The Galerkin Finite Element Method is used to transform 
the fourth-order differential equation governing static 
displacement of the beam into a system of linear algebraic 
equations. In order to maintain continuity of slopes and 
displacements from element to element, a family of cubic shape 
functions is introduced [{ref. 5 and ref. 6]. An approximate 
solution for lateral displacement is defined as: 

vs § ="? 
where: 


v = exact solution in continuous space (23) 
¥ = approximate solution in discrete space 


= 


= vector containing shape functions 
= vector containing lateral displacements and slopes 


<u 


The next step is to form a measure of the error of the 
approximation, or residual: 
& salgi?') - p,(s) (24) 
where p, is the lateral load and s is the independent 


variable. Substituting equation (23) into equation (24) 


yields the following equation: 
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R= [EI(N7H) "|" - py(s) (25) 


The residual can be minimized if: 
[ fiR(s) ds = 6 (26) 
D 


where the quantity on the right hand side of the equation is 
the null vector. Substitution of equation (25) into equation 


(26) gives the following result: 
[ Nez a9) "ds - [ Np, (s) ds = 0 (27) 
D D 


Integration by parts is performed twice on equation (27), 


resulting in: 


N[EI(N7#)")'|, - (MEI (N77) "|, 
xo ff 37-2. iy 7 "al 
+ fn EI(N’V) "ds - | Hp, (s) ds = 0 


where |, indicates an evaluation at the boundary points of the 
structure. Since the displacement vector is a constant, 


equation (28) can be rewritten as: 


N(EI(N")"'3|, - (M)/EI(N?) "FI, 
a I Ty if y; 6 ae 
+ fun EI(N*) "dst - [ fie, (s) ds = 6 


From the beam equilibrium equation, the shear force, V, is 
defined as: 
V = EIvll (30) 


and the moment, M, by: 
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M = EIv" (31) 
Therefore, the boundary term load vectors are defined as: 
V = N[EI(N*)")'?|, (32) 
and: 
M = (N)'EI(N*)"F|,, (33) 


A system stiffness matrix is defined as: 


K; = f (N) “EI (N7) "dg (34) 
D 
and a system force vector by: 
F, = f{ Np,(s) ds (35) 


Substitution of equations (32) through (35) into equation (29) 
yields the following system of linear equations: 


= = 


Vi, - Ml, + K,0- F, = 0 (36) 
If a load vector of internal and external applied loads is 
defined as: 
F, =F, + Ml, - A, (37) 
equation (36) reduces to: 


K,¥ = F, (38) 
The global bending stiffness matrix is constructed from the 
union of the elemental bending stiffness matrices and the 
global bending force vector is constructed from the union of 


the elemental bending force vectors. 


Zo 


C. FINITE ELEMENT BAR EQUATION DEVELOPMENT 

The development of the Galerkin Finite Element Method for 
the bar equilibrium equation is done in a manner similar to 
that for the beam equation. However, since the bar equation 
is a first-order differential equation, a family of linear 
shape functions can be used to maintain continuity of 
displacement from element to element [ref. 5 and ref. 6]. 
Once again, an approximate solution for the axial 


displacement, u, is formed as follows: 


us i = NG 
where: 
u = exact solution in continuous space (39) 
& = approximate solution in discrete space 
N = vector containing shape functions 
i = vector containing axial displacements 
The residual is defined as: 
R(s) = [AE(N7@)']! + p.(s) (40) 
The residual can be minimized if: 
[ NIR(s)1ds = 6 (41) 
D 


Substituting equation (40) into equation (41) results in the 


following equation: 
[ fitaz ia) ids + f Np,(s) de = 0 (42) 


Unlike the beam equation development, only one integration by 


parts is performed to yield: 
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AEN (NT)'|, - [ah taza) ds + [ fp. (s) ds =G (43) 


where }, indicates an evaluation at the boundary points of the 
structure. Since the axial displacement vector is a constant, 


equation (43) can be rewritten as: 
N( AEN) ‘al, - [a TAB) dst + | Nip.ds = 9 (44) 


From the bar equilibrium equation, the axial force, F, is 


defined as: 


F = AEu!' (45) 


Therefore, the boundary term load vector vectors are defined 


as; 


B = AEN(H) dl, (46) 


A stiffness matrix is defined as: 
K, =f (M)'(AE(N")] ds (47) 
and a system force vector by: 
F, = [ Np.(s) ds (48) 


Substituting equations (46) through (48) into equation (44) 


yields the following equation: 


B-K,a+ F,=0 (49) 


By defining a new load vector as: 


Z5 


Fo = Fer (50) 


equation (49) reduces to: 


K,a = F, (51) 

The global axial stiffness matrix is constructed from the 
union of the elemental axial stiffness matrices and the global 
axial force vector is constructed from the union of the 


elemental axial force vectors. 


D. THE ELEMENTAL STIFFNESS MATRIX 

The global Galerkin FEM stiffness matrices in equations 
(38) and (51) are constructed from the union of elemental 
axial and bending stiffness matrices. A single beam element 
has four degrees of freedom, which yields a 4x4 elemental 
stiffness matrix. A single bar element has two degrees of 
freedom, which yields a 2x2 elemental stiffness matrix. These 
elements can be combined into a single bar-beam element (fig. 


4.1) which has six degrees of freedom [ref. 5 and ref. 6]. 





Figure 4.1: Bar-Beam Element - Degrees of Freedom 
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This results in a 6x6 elemental stiffness matrix. The 


elemental displacement vector can be expressed as: 


GA? = <0," 62 98,7 09") 8,7',/8,2°> 
where for the i“ element: 


8,1’ = axial displacement at local node 1 

§,/' = lateral displacement at local node 1 

8,/’ = beam slope at local node 1 oe 
8,4 = axial displacement at local node 2 

$,?’ = lateral displacement at local node 2 


8,2’ = beam slope at local node 2 


The elemental force vector can be expressed as: 


ee ie re ee cs ee t’> 
where for the i™ element: 


i/ 
f, 


axial force at local node 1 


i/ 
f, 


lateral force at local node 1 
(53) 
moment at local node 1 


i/ 
f, 


f,/' = axial force at local node 2 


lateral force at local node 2 


i! 
f, 


f,/’ = moment at local node 2 


Therefore, the combination of the Galerkin bar and beam 


equations yields: 
k2'§ 4! = fa! (54) 


where the elemental stiffness matrix takes the following form: 
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AE _ AE 
7 0 0 TL 0 0 
EI . EI EI . EI 
0 22 2 eee 6 ZL 
i~p i? ip 1; 
0 6f. “eo Wc ieee 
=! 1; fy wy 1; 
si ad Wa: AE oy 
~ se 6 2 0 0 
Jy 1; 
ye he EI EI 
Cum aie oe Ty 2fi _,#L 
Family iy gullies 
EI EI er ee 
eG en ac q=t 
Dae 1, ae by 


It is apparent from the form of the elemental stiffness matrix 


that the bar and beam have uncoupled behavior. 


E. COORDINATE TRANSFORMATION OF THE ELEMENTAL SYSTEM OF 
EQUATIONS 

Each element in the arch will have a unique orientation 
with respect to the global x and y axes. In order to solve 
the global system of equations, it is necessary to transform 
each elemental Galerkin equation into global coordinates. The 
global reference coordinate system is defined as_ the 
horizontal and vertical axes of the entire arch. Each element 
makes an angle a, with the horizontal axis and an angle &; with 
the vertical axis. 

The local displacements and forces (indicated by prime 


superscript) are defined as follows: 
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6,2" = §,7cos(a,) + 5,7cos(B,;) 


6,7’ = -8,4cos(B,) + 6,4cos(a@,) 
6.7 = 5.1 
3 3 (56) 
8,2" = 8,4cos(a;) + 6,4cos(B,) 
6,7’ = -6,/cos(B,) + &,%cos(a,) 
7° = 8,7 
and: 
£7 =e6 cosi( Ga) ee £2c0s (p,) 
£,"" = -f,!cos8(B,) + f,4cos(e,) 
f.i' = Fi 
3 3 (57) 
ft’ = f,4cos(a,) + f,4cos(B,;) 
f,/' = -f,4cos(B;) + f57cos(a,) 
ce ss f,4 
Therefore a transformation matrix can be defined as: 
cos({a,) cos(B,;) 0 0 0 0 
-cos(B,;) cos(a@,;) 0 0 0 0 
T « 0 0 iL 0 0 0 (58) 
0 0 0 cos({a;) cos(B,;) 0 
0 0 0 -cos(B,) cos(a,) 0 
0 0 0 0 0 a 
This reduces equations (56) and (57) to: 
bi! = Figs (59) 
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and: 


£ al462 (60) 

where: 
(87)? = ¢8,7,857,6,,/,057,8.7,8,7> (61) 
(E*)? = <6, 2s? ee ee (62) 


The transformed elemental stiffness equation becomes: 


x Di§S= Pit (63) 
Multiplying both sides of equation (63) by the inverse of the 


transformation matrix, which is an orthogonal matrix, yields: 


(4) 72 (24) ae? (64) 
where the elemental stiffness matrix in terms of global 


coordinates is defined as: 


ki = (PD?) tk?’ (TP?) (65) 





Figure 4.2: Coordinate Transformation 
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F. SOLUTION 

The global stiffness matrix and global force vector are 
constructed from the union of the elemental stiffness matrices 
and elemental force vectors, which yields the following 
equation: 


KA =F 
where: 


K = global stiffness matrix (66) 

A = global displacement vector 

F = global force vector 

By multiplying both sides of the equation by the inverse 

of the global stiffness matrix, the global displacement vector 
can be calculated. The global displacements are then 
transformed back into the local displacements and can be used 
in conjunction with the elemental stiffness matrix to 


calculate the local forces: 
k2'§ 2! = FA’ (67) 
These local forces and bending moments are used to 
calculate the stress at each nodal point. The stresses of 
internal global nodal points are averaged since local nodal 


point 2 of the it" element is the same as local nodal point 1 


of the (i+1)*® element. 
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V. PROGRAM DESCRIPTION 
The Fortran 77 code originally used in references 5 and 6 
to analyze circular arches was modified to analyze the non- 
circular arches of interest in this study. A copy of the 
program is included in Appendix B. The program reads 
information from an input file which provides the following 


user-supplied data: 


L = horizontal distance spanned by arch 

H = vertical distance spanned by arch 

HGT = depth of cross-section 

YOUNG = Young’s Modulus 

YIELD = yield strength 

NEL = number of elements 

ISTRAT = ADS parameter which designates strategy 

IOPT = ADS parameter which designates optimizer 

IONED = ADS parameter which designates one-d search method 
IPRINT = ADS print control parameter 


IGRAD = ADS parameter which designates method of gradient 
calculation 


DV_BG = initial value of design variable 


DV_LO = lower side constraint on design variable 


upper side constraint on design variable 


DV_UP 
CLAN = node at which concentrated load is applied 
FX = concentrated horizontal force (positive to the right) 


in pounds 
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FY = concentrated vertical force (positive upwards) in 
pounds 


FM = concentrated moment (positive counter-clockwise) in 
pound-inches 


FA = distributed load in pounds per inch 


BX_ = boundary condition on horizontal displacement at 1, 
first node, or 2, last node (= 1 - fixed,= 0 - free) 


BY = boundary condition on vertical displacement at l, 
first node, or 2, last node (= 1 - fixed,= 0 - free) 


BM_ = boundary condition on rotation at 1, first node, or 2, 
last node (= 1, - fixed,= 0 - free) 


LABEL = string character identifying particular case study 


The program then utilizes several subroutines to perform the 
FEM analysis and call the ADS program to perform the 
optimization. 

An outline of the program, ARCHOPT, and its subroutines is 
provided in fig. 5.1. The first subroutine called by the main 
program is OPTIMIZATION TOOL, which establishes the ADS 
parameters before the first call of the ADS program. The 
first call serves to override some of the default parameters 
in order to "fine-tune" the program. After ADS is called, 
OPTIMIZATION TOOL calls subroutine EVAL to evaluate the 
objective function and constraints, which are functions of the 
design variables. 

Subroutine EVAL calls subroutine ARCH STRESS, which in 
turn calls subroutines FORM and FORCE _VECTOR. The latter 


subroutines are used to form the global stiffness matrix and 
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the global force vector, respectively. Subroutine BNDRY 
imposes the user-supplied boundary conditions by modifying the 
stiffness matrix and the force vector. An equation solver 
from the IMSL library, L2ARG, is used to solve the Galerkin 
FEM equation by matrix inversion. Once the_ global 
displacements are known,they can be used to calculate the 
nodal stresses, which is done in subroutine STRESS. Once the 
design constraints are evaluated for the initial design, the 
problem is returned to ADS, where an updated design is chosen. 
This process continues until the solution converges. 

Once the termination criteria for optimization is 
satisfied, the main program calls subroutine ARCH _OUT, which 
creates an output file, ARCHOUT, which contains the problem 
parameters, the optimized design, and the arch volume. Copies 
of the output files for all case studies are provided in 
Appendix D. 

The methodology of the program was thoroughly validated in 
references 5 and 6. Once the modifications were made, the 
program was used to analyze several statically determinate 
structures (i.e. a cantilever beam) to verify that the 
approximate values of stress and displacement came close to 
the exact values. The results confirmed that the FEM analysis 
provided an excellent approximation, with less than 0.8% error 
for a twelve element model of a cantilever beam with a 


lateral concentrated load at the free end. 
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Program 
ARCHOPT 


OPTIMIZATION 
FOL 


File 
ARCHIN 


ARCH our 





Figure 5.1: Program ARCHOPT Structure 
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VI. CASE STUDIES 

The optimized designs for eleven different cases of 
loading and boundary conditions are presented. Cases (1) 
through (4) are hinged-hinged arches which span equal 
distances in the horizontal and vertical directions. Each 
case has a different loading. Cases (5) through (8) are 
arches which span equal horizontal and vertical distances, but 
which have different boundary conditions. Cases (9) and (10) 
are hinged-hinged arches which span different distances in the 
vertical direction. Plots of arch shape, base dimensions, and 
stress are provided for each case. It should be note that the 
distributed loads are applied in a direction normal to each 
element. This would closely approximate a situation where the 
arch was under hydrostatic, or pressure loading. Twelve 
elements were used to model the arch in all cases. This 
implies that there are thirteen nodal points and thirteen 
stress constraints. A Young’s Modulus of 30,000,000 psi and 
a S, of 52,000 psi are used in all cases. The output files 


from program ARCHOPT are presented in Appendix C. 
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A. CASE 1: HINGED-HINGED ARCH WITH DISTRIBUTED LOAD 

This design has nine active stress constraints and 
thirteen active side constraints (lower limit on base). The 
stress due to axial loading is higher than the bending stress 
at all nodes. This indicates that the minimum area cross- 
section is better able to withstand an axial load than a 
bending moment. A predominantly axial stress distribution is 
more efficient since the entire cross-section is stressed to 
the same degree. In the case of a bending stress 
distribution, the material has to be strong enough to 
withstand the extreme fiber stress, which means that the 
internal part of the cross-section will be stressed below the 


yield strength. 


@e H = 32.0 inches 


32.0 inches 


eL 


eh = 1.5 inches 


@ Distributed load = -100 pounds/inch 
@ Volume = 10.324 cubic inches 
loo /q 


Figure 6.1: Case 1 - Initial Shape 
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Case 1: Arch Shape and Width 
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Figure 6.2: Case 1 - Arch Shape and Width 
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Case 1: Stress 
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Figure 6.3: Case 1 - Stress 
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B. CASE 1A: HINGED-HINGED BEAM WITH DISTRIBUTED LOAD 

This case is presented for comparison to Case 1. The 
slopes of the elements are not allowed to vary, so the arch 
remains a straight beam. There are seven active stress 
constraints and no active side constraints. The bending 
stresses dominate and the axial stresses are negligible. The 
volume of this arch is more than six times greater than that 


for Case l. 


e@ H = 32.0 inches 

e L = 32.0 inches 

eh = 1.5 inches 

@ Distributed Load = -100 pounds/inch 


@ Volume = 68.387 cubic inches 


j60 Vb) 


Figure 6.4: Case 1A - Initial Shape 
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Case 1A: Arch Shape and Width 
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Figure 6.5: Case 1A - Arch Shape and Width 
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Case 1A: Stress 
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Figure 6.6: Case 1A - Stress 
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C. CASE 2: HINGED-HINGED ARCH WITH CONCENTRATED LOAD 

This arch has a concentrated vertical load applied 
downward at node 7. There are seven active stress constraints 
and thirteen active side constraints (lower limit on base). 
The arch has taken the shape of two straight lines meeting at 
the center node. At eight of the thirteen nodes the axial 


stress dominates. 


e H = 32.0 inches 
e L = 32.0 inches 
eh = 1.5 inches 


@ Concentrated Load = 2000 pounds (downward) 


@e Volume = 10.239 cubic inches 


2600 Ib 


Figure 6.7: Case 2 - Initial Shape 
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Case 2: Arch Shape and Width 
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Figure 6.8: Case 2 - Arch Shape and Width 
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Case 2: Stress 


-o, 
oom 
”) 
joe 
al 
2] 
”—) 
® 
Sug 
cline? 
ep) 





Figure 6.9: Case 2 - Stress 
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D. CASE 3: HINGED-HINGED ARCH WITH CONCENTRATED LOAD 

This design has five active stress constraints and twelve 
active side constraints (lower limit on base). The arch has 
taken the shape of two straight lines meeting at node 7, where 
the concentrated load is applied. Once again, this design 


shows a preference for axial loading over bending moments. 
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Figure 6.10: Case 3 - Initial Shape 
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Case 3: Arch Shape and Width 
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Figure 6.11: Case 3 - Arch Shape and Width 
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Case 3: Stress 
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Figure 6.12: Case 3 - Stress 
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E. CASE 4: HINGED-HINGED ARCH WITH CONCENTRATED MOMENT 
This design has one active stress constraint and twelve 
active side constraints (lower limit on base). The arch 
assumes the shape of a straight line, which is the minimum 
shape for a feasible design. The maximum base is at node 7, 


where the concentrated moment is applied. 


@e H = 32.0 inches 


32.0 inches 
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Figure 6.13: Case 4 - Initial Shape 
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Case 4: Arch Shape and Width 
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Figure 6.14: Case 4 - Arch Shape and Width 
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Case 4: Stress 
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Figure 6.15: Case 4 - Stress 
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F. CASE 5: ROLLER-HINGED ARCH WITH DISTRIBUTED LOAD 

This design has ten active stress constraints and twelve 
active side constraints (lower limit on base). The arch has 
assumed a nearly circular shape. In this case, the bending 
stress dominates throughout the arch. It appears that the 
arch has attained the maximum degree of axial loading that can 


be achieved with these boundary conditions. 


e H = 32.0 inches 


e L = 32.0 inches 

eh = 1.5 inches 

@ Distributed Load = -100 pounds/inch 
@ Volume = 11.329 cubic inches 


ISO lin 


Figure 6.16: Case 5 - Initial Shape 


52 


Case 0: Arch Shape and Width 
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Figure 6.17: Case 5 - Arch Shape and Width 


oS 


Case 5: Stress 
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Figure 6.18: Case 5 - Stress 
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G. CASE 6: ROLLER-FIXED ARCH WITH DISTRIBUTED LOAD 
This design has four active stress constraints and ten 
active side constraints (lower limit on base). The bending 


stresses dominated for the most part, but the axial stresses 


are significant. The maximum base is at the fixed end, as 
expected. 
e H = 32.0 inches 


eL 32.0 inches 


eh = 1.5 inches 
@e Distributed Load = -100 pounds/inch 


@ Volume = 14.431 cubic inches 


106 \oi;, 


Figure 6.19: Case 6 - Initial Shape 


55 


Case 6: Arch Shape and Width 
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Figure 6.20: Case 6 - Arch Shape and Width 
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Case 6: Stress 
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Figure 6.21: Case 6 - Stress 
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H. CASE 7: FIXED-FREE ARCH WITH DISTRIBUTED LOAD 

Using the cross-sectional depth of the previous case 
studies resulted in a design that consistently violated the 
stress constraints at the nodes at and close to the fixed end 
of the arch. Therefore a cross-sectional depth of 3 inches 
was used. The resulting design has seven active stress 
constraints and six active side constraints (lower limit on 
base). The arch has taken the shape of a tapered beam with a 
maximum base at the fixed end. The bending stresses dominate 
and the axial loads are negligible. Because of the lack of 
support at the free end, the arch is unable to attain axial 


loading in the elements. 


eu 32.0 inches 


eL 32.0 inches 


eh = 3.0 inches 
@ Distributed Load = -100 pounds/inch 


@® Volume = 72.614 cubic inches 


}O9 ia 


Figure 6.22: Case 7 - Initial Shape 
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Case 7: Arch Shape and Width 
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Figure 6.23: Case 7 - Arch Shape and Width 
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Case 7: Stress 
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Figure 6.24: Case 7 - Stress 
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I. CASE 8: FIXED-FIXED ARCH WITH DISTRIBUTED LOAD 

This design has six active stress constraints and eleven 
active side constraints. The arch curves outward slightly in 
the center. With the exception of the nodes at the ends, 
axial stresses dominate in the arch. The maximum base 


dimensions occur at the fixed ends. 


e H = 32.0 inches 

e L = 32.0 inches 

eh = 1.5 inches 

e Distributed Load = -100 pounds/inch 


@e Volume = 11.308 inches 


100 Wesin 


Figure 6.25: Case 8 - Initial Shape 
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Case 8: Arch Shape and Width 
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Figure 6.26: Case 8 - Arch Shape and Width 
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Case 8: Stress 
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Figure 6.27: Case 8 - Stress 
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J. CASE 9: HINGED-HINGED ARCH WITH DISTRIBUTED LOAD 

This case differs from Case 1 only in the distance spanned 
in the vertical direction , H. In this case, H = 18.475 
inches. There are eleven active stress constraints and 
thirteen active side constraints (lower limit on base). The 
arch has assumed a curved shape similar to Case l. Once 


again, there is a clear preference for axial loading. 


e H = 18.475 inches 

e L = 32.0 inches 

eh = 1.5 inches 

@ Distributed Load = -100 pounds/inch 


@ Volume = 8.391 cubic inches 





Figure 6.28: Case 9 - Initial Shape 
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Case 9: Arch Shape and Width 
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Figure 6.29: Case 9 - Arch Shape and Width 
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Case 9: Stress 
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Figure 6.30: Case 9 - Stress 
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K. CASE 10: HINGED-HINGED ARCH WITH DISTRIBUTED LOAD 

This case differs from the previous case only in the 
distance spanned in the vertical direction, H. There are 
eight active stress constraints and twelve active side 
constraints. Once again, the arch curves outward and the 


axial stresses dominate. 


@ H = 55.426 inches 

e L = 32.0 inches 

eh = 1.5 inches 

@ Distributed Load = -100 pounds/inch 


@® Volume = 15.149 cubic inches 





Figure 6.31: Case 10 - Initial Shape 
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Case 10: Arch Shape and Width 
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Figure 6.32: Case 10 - Arch Shape and Width 
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Case 10: Stress 
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Figure 6.33: Case 10 - Stress 
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The 


VII. CONCLUSIONS 


conclusions of this study are: 


The objective function appears to have local optimum 
points, because the initial shape of the arch affected the 
output from ADS. This required a trial-and-error approach 
to find the true optimal design for each case. 


In the optimum designs, the base dimensions were reduced 
to the maximum extent possible, and the arch took a 
rounded shape to decrease bending moments and increase 
axial loading. One exception to this was the case of the 
fixed-free arch, which was unable to attain axial loading 
and took the shape of a straight, tapered beam. The only 
other exception was Case (4), the hinged-hinged arch with 
a concentrated moment. 


The smooth, curved shapes that most of the arches assumed 
confirmed that modelling an arch as a series of straight 
elements was a valid approach as long as the number of 
elements was relatively high (greater than eight). 


When using the combination of strategy and optimizer 
(S.L.P. and M.M.F.D.) that this investigation employed, it 
is critical that the initial design be feasible. Since 
the initial shape in all cases was a straight beam, the 
initial base dimensions had to be large enough to allow 
the beam to be able to withstand the bending stresses that 
would dominate in such a structure. If the initial design 
was not feasible, the resulting output from ADS had many 
violated constraints. 


The most inefficient method of spanning the distance from 
point A to point B was the cantilever beam (Case (7)). 


Case Studies (5) through (8) confirmed the conclusions of 
references 5 and 6 that, for a given loading, the more 
statically indeterminate the structure, the greater the 
efficiency. However, Case (1) contradicted this 
conclusion since it is less statically indeterminate than 
Case (8) but has a slightly smaller volume. 
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Suggested areas for future research are: 


e Optimization of an arch with a more complex cross-section 
(i.e. thin-walled tubes, WF’s, etc.) 


e Optimization of arches with geometric constraints which 
prevent global buckling and local crippling. 
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APPENDIX A: JUSTIFICATION FOR OMITTING SHEAR STRESSES 


(The following appendix is taken from ref. 5) 


The shear stress distribution through a beam of 


rectangular cross-section has a parabolic distribution along 


the depth of the member. The maximum shear stress, located at 


the neutral axis of the beam, is (from ref. 10): 


where: 


€ = maximum shear stress 


mar 
V = shear force 
A = cross-sectional area 


(B.1) 


The maximum normal stress due to bending is given by: 


where: 
0, = maximum normal stress 


M = bending moment 

c = distance from N.A. to extreme fiber 
3 

I= a = moment of inertia 


Redefining the normal stress in terms of the 


sectional dimensions yields: 
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(B.2) 


cross~- 


0, = M(h/2)/(bh?/12) 
or: (Be3') 
0, = 6M/hA 


The ratio of the maximum shear stress to the maximum 


normal stress due to bending is denoted by r and given by: 


— — (B.4) 





Substituting equations (B.1) and (B.3) into equation (B.4) 


yields: 


r= (1.5V/A) /(6M/BA) 


Vh 


Tri=-—— 
4M 


For the cases investigated in this study, the maximum value r 
can attain is when the loading is that of a uniformly 


distributed load, p,. In this case: 


V = pL 


and: (B.6) 
M = p,L?/2 


which upon substitution into equation (B.5) yields: 
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r= (p,L) h/4(p,L?/2) 


prt (B.7) 
Be 
2L 


I= 


The use of the beam equation requires that length of the 
beam to be at a minimum ten times the height, or in other 


words: 
L210h (B.8) 


To maximize the value of r, let L equal 10h, the minimum 
allowable length. Substituting this value of L into Equation 
(B.7) yields: 

resehf2(10h) 

or: (B.9) 

Tels 20 
Therefore, the maximum shear stress accounts for less than 

5% of the bending stress developed in the structure. Five 
percent is high considering this analysis over-assumed the 
value of the shear stress by assigning the maximum shear 
stress to the entire cross-section of the beam. Moreover, at 
the outermost fibers where the normal stress is a maximum, the 
shear stress is zero. Therefore, under the circumstances of 
this study, the addition of shear stresses was deemed to be 


unwarranted. 
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APPENDIX B: Program ARCHOPT 
FILE: ARCHOPT FORTRAN A 


PROGRAM ARCHOPT 
(C36 3 36 3 36 3 3 36 3 2 3 3 2 3 3 2 3 38 2 2 2 38 2 2 D6 2 2 2 3 222 3 2 2 2 36 2 2 2 28 2 3 30 2 2 2 3 2 36 2 2 38 36 2 2 DE 2 2 3 3 


Cx x 
Cx ARCH OPTIMIZATION ANALYSIS CODE x 
Cx x 


C36 Be BE BE FE IE UE HE DE HE EE EE DEE DE DE DE SE DE SE OE DE SE SE SE SE SE DE EE SE OE DE DE DEE DE DE SE SE DE OE SE DE SE DE OE OE OE DE DE SE OE DEE OE DE DE DE DE DE DE DE IE IE IE XE 


ALPHA.... TRANSFORMATION ANGLE OF ELEMENT (ANGLE TO X-AXIS) 
BAVE.....-THE AVERAGE BASE DIMENSION ACROSS AN ELEMENT 
BASE....+ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS 
BASEL...-ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS LOWER 
SIDE CONSTRAINT 
BASEU....ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS UPPER 
SIDE CONSTRAINT 
BETA ....TRANSFORMATION ANGLE OF ELEMENT (ANGLE TO Y-AXIS) 
B_1.....-BOUNDARY TERMS APPLIED AT END "1" 
B_2......BOUNDARY TERMS APPLIED AT END "2" 
Cl1,..,CS5.CONSTANTS RELATED TO ELEMENT STIFFNESS COEFFICIENTS 
CLAN....-CONCENTRATED LOAD APPLICATION NODE (THE NODE FX,FY,FM ARE 
APPLIED 
DOF......DEGREE OF FREEDOMS (UNKNOWN DISPLACEMENTS & SLOPES) 
DVIBG....DESIGN VARIABLE #1 (BASE DIMENSION) INITIAL ESTIMATE 
DVILO....DESIGN VARIABLE #1 (BASE DIMENSION) LOWER SIDE CONSTRAINT 
DV1UP....DESIGN VARIABLE #1 (BASE DIMENSION) UPPER SIDE CONSTRAINT 
DV2EBG....DESIGN VARIABLE #2 (SLOPE) INITIAL ESTIMATE 
DV2LO....DESIGN VARIABLE #2 (SLOPE) LOWER CONSTRAINT 
DV2UP....DESIGN VARIABLE #2 (SLOPE) UPPER CONSTRAINT 
EK...++--6X6 ELEMENT STIFFNESS MATRIX IN LOCAL X,Y COORDINATES 
EKPR.....6X6 ELEMENT STIFFNESS MATRIX IN ELEMENT LOCAL COORDINATES 
ELEN.....LENGTH OF ELEMENT 
Hsls ©1610. % 10 FORCE VECTOR OF SYSTEM 
FA.....--CONSTANT DISTRIBUTED LOAD OUTWARD FROM END TO END 
FM....2ee CONCENTRATED MOMENT AT CLAN 
FX.....-+eCONCENTRATED LOAD IN X DIRECTION AT CLAN 
FY......-CONCENTRATED LOAD IN Y DIRECTION AT CLAN 
G.eeeeeee THE ARRAY OF CONSTRAINT FUNCTIONS 
GAMMA....6X6 ELEMENT TRANSFORMATION MATRIX 
+ ayiane “ois - (NDOF JX(NDOF) GLOBAL STIFFNESS MATRIX 
HGT.....-CONSTANT DEPTH OF CROSS-SECTION 
INFO....-ADS PARAMETER USED TO SIGNAL THAT THE OPT IS COMPLETE 
IPRINT...ADS PARAMETER USED SELECT THE DATA OUTPUT FORMAT 
ITERATE..THE NUMBER OF TIMES ADS IS TO BE RELOADED WITH THE 
PRECEEDING DATA 
IWK. 2.20 ADS INTERNAL WORK SPACE ARRAY 
NCON.....-NUMBER OF DESIGN CONSTRAINTS 
NDOF....-NUMBER OF DEGREES OF FREEDOM 
NDV.....-NUMBER OF DESIGN VARIABLES 
NEL.....-NUMBER OF ELEMENTS 
NRIWK...-ADS INTERNAL WORK SPACE ARRAY DIMENSION 
NRWK....eADS INTERNAL WORK SPACE ARRAY DIMENSION 
NSNP.....NUMBER OF SYSTEM NODAL POINTS 
OBJ......-THE OBJECTIVE FUNCTION OF THE OPTIMIZATION 
OPTDCS...OPTIMIZATION DECISION TO OPTIMIZE THE PROBLEM OR NOT 
P1...P5..PARAMETER DIMENSION CORRESPONDING TO THE NEL» NSNP, NCON, 
NDOF, AND NDV RESPECTIVELY 
PRCSN....-THE PRECISION DESIRED TO SOLVE THE FEM SYSTEM OF EQUATIONS 
SIGMA_B..THE ELEMENTAL NORMAL STRESS DUE TO BENDING 
SIGMA_N..THE ELEMENTAL NORMAL STRESS DUE TO AXIAL FORCES 
SIGMA_T..THE MAXIMUM TOTAL STRESS IN EACH ELEMENT 
SX.eee-eeeGLOBAL HORIZONTAL COORDINATE 
SY..-ee-eGLOBAL VERTICAL COORDINATE 
U.seeeeeee THE "DISPLACEMENT" VECTOR OF THE SYSTEM OF LINEAR EQUATIONS 
VLB...+++ADS ARRAY CONTAINING UPPER SIDE CONSTRAINTS 
VUB....--ADS ARRAY CONTAINING LOWER SIDE CONSTRAINTS 
WK..-e+2+-eADS INTERNAL WORK AREA 
Xe eeeeeeeADS ARRAY CONTAINING THE VALUES OF THE DESIGN VARIABLES 
YIELD....VYIELD STRENGTH OF THE ARCH MATERIAL 
YOUNG....YOUNG'S MODULUS OF THE ARCH MATERIAL 
C36 3 DE 3 ES IE 3G Se DE DE DE DE DE DE OE DE DE DE DE OE DE DE DE OE DE DE DE DE DE DE DE DE 3 DE 3 DE OE DE DE DE DE DE DE DDE 3 DE DE DE DE DE 3 DE DE DE DE 3 DE DE DE OE DF DE DE DE DE DE DE OE 
C234567 eee + DECLARE THE NUE MES SIMeT sc. 4 ole) cc's ee felace's 6 6.0 6. 6.6.6 6 60 0 6 0 0 0 Shenetee 46 
INCLUDE 'ARCHCOM FORTRAN’ 


ADQOOAOAOOIOIAADAAAAIAAVOAVANIVDAAAIAAVAANAIAIVDANIAIDIAIAIAINIDIAIGDININIAIANAARANANAIAINANIANNAAANAAQD|AN0NNO 
o 
“A 


15) 


FILE: ARCHOPT FORTRAN A 


Ec read the input Parameters. .. .oosccsc ccs ss eleenetes clmumenINCDS |. 
OPEN(8, FILE='ARCHIN', STATUS='OLD') 


READ(8,*) L,»,H,HGT 

READ(8,*) YOUNG, YIELD 

READ(8,*) NEL,»,ISTRAT, IOPT, TIONED, IPRINT,IGRAD 
READ(8,x*%) DV1BG,DV1L0O,DV1UP 

READ(8,%) DV2BG,DV2L0,DV2UP 

READ(8,x%) CLAN,FX,FY 

READ(8,x*%) FM,FA,OPTDCS 

READ(8,*) ITERATE,PRCSN,BX1,BY1,BM1 
READ(8,x*%) BX2,BY2,BM2 

READ(8,x%) LABEL 


--- define CONSTANTS. oS cbc ce baleen 0 ote cleletete so ccs slseetcreiete cise ss oo suai 
NSNP = NEL + 1 
NDOF = 3*¥NSNP 
NCON = NSNP 
NDV = NSNP + (NEL = 1) 


a0 


«+e OPTIMIZE THE PROB LEMS ..27e < « oie 0 eles <0 0 she 0.0: sleie «clele ecctereterenn > « clams 
CALL OPTIMIZATION_TOOL 


eee e- COMPILE AND FORMAT THE GCUTPURG s.. os « SRNR as oS I ees sé cut 
CALL ARCH_OUTPUT 


END 
HE HE HE HE HE HE HE IESE HE HE IE SE HE HE SE SE SE HE SE SE HE BE SE HE HE SE SE HE HE SE IE HE EE HE HE HE HE EE EK HE HE HE HE HH EK RH MK KR KKK RR KKM ¥ 


00 09 00a 200 


SUEROUTINE OPTIMIZATION_TOOL 


THIS SUBROUTINE DIRECTS THE PROGRAM FLOW OPTIMIZATION DECISION 
I.E.» OPTIMIZE THE PROBLEM OR NOT. IT ALSO SERVES TO SET UP & 
EXECUTE THE ADS OPTIMIZATION SOFTWARE. 


«ee Gdeclare the VAPLEDLES. 0. eee eee se ete csscerttsercececceoess 
INCLUDE ‘ARCHCOM FORTRAN’ 
INTEGER I 
DO 100 I=1,NSNP 

BASE(I) = DV1IBG 
BASEL(TI) DV1LO 
BASEU(T ) DV1UP 
100 CONTINUE 
DO 150 I=1,NEL-1 
SLP(I) = DV2BG 
SLPL(I) = DV2LO 
SLPU(I) = DV2UP 
150 CONTINUE 
SLP(NEL)=DV2BG 
SLPL(NEL)=DV2LO 
SLPU(NEL)=DV2UP 


Cc ««.eCOMBINE BASE AND SLP ARRAYS INTO DESIGN ARRAY......... 
DO 200 I=1,NSNP 
X(I) = BASE(T) 
VLB(IT) = BASEL(T) 
VUB(I) = BASEU(T) 
200 CONTINUE 
DO 250 J=NSNP+1,NDV 
X(J) = SLP(J-NSNP ) 
VLB(J)=SLPL(J-NSNP ) 
VUB(J)=SLPU(J-NSNP ) 
250 CONTINUE 
Cc 


c eee MAKE OPTIMIZATION DECISION sisiecs « is ss ole ororeteieiets o1cccleis 6c <Ber. 
IF (OPTDCS .NE. 1) THEN 
CALL EVAL 
RETURN 
ENDIF 


o+++DEFINE THE SIZE OF THE WORK ARRAYS FOR ADS.....ceceseccees 


Qgaqgaqgng 00 


a0 
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FILE: ARCHOPT FORTRAN A 


NRA=65 
NCOLA=97 
NRWK = 40000 
NRIWK = 2000 
oo ready to SemeIINIMIE EMEC. fc slclsle cislsleise esse ccc cece eee eee eee ecess 
DO 280 I=1,NSNP 
IDG(I)=0 
280 CONTINUE 
INFO =-2 
_CALL ADS CINFO,ISTRAT, IOPT,IONED, I PRINT, IGRAD,»NDV,NCON,X,VLB,VUB, 
es OBJ,G,IDG,NGT,IC, DF,A,NRA, NCOLA, WK» NRWK; IWK, NRIWK) 
“IMKC = 


ae CALL ADS CINFO,ISTRAT,IOPT,IONED, ITPRINT, IGRAD,»,NDV,NCON, X,VLB,VUB,> 
OBJ,G,IDG,NGT,IC,DF,A,NRA, NCOLA,WK, NRWK, IWK, NRIWK) 


ok 


~+-e-eValuate the objective function and constraints.........6. 
IF (INFO .NE. 0) THEN 
CALL EVAL 


C3 BE BE IE I DE DE IE HE HE IE IE IE IE IE DE DE DE DE IE DE IE HE DE IE DE DE IE IE IE IE IE IE IE IE IE IE IE IE IE IE IE DE IE IE IE IE DE IE IE IE IE DE IE IE DE DE IE IE DE IE DE IE DE DE IE IE DE IK 


SUBROUTINE EVAL 


C THIS SUBROUTINE IS USED TO EVALUATE THE OBJECTIVE FUNCTION, 
C Seer tee e FUNCTIONS, AND SIDE CONSTRAINTS OF THE OPTIMIZATION 
c OBLENM. 


c Teweaecrare the Variables... cceesceses Reunite! 6 < SreTkteo es ss aa etene 
INCLUDE "ARCHCOM FORTRAN’ 
INTEGER I,J>K 


REAL PI 
" PARAMETER (PI=3.141593) 
Cc eee +SEPARATE THE DESIGN ARRAY INTO BASE AND SLP ARRAYS 


DO 50 I=1,NSNP 
BASE(I) = X(T) 
50 CONTINUE 
DO 75 I=1,NDV-NSNP 
SLP(I) = XC I+NSNP) 
75 CONTINUE 
DX=L/NEL 
SY(1)=0.0 
DO 80 I=1,NEL-1 
SYCI+1)=DX¥SLP(I) + SY(T) 
80 CONTINUE 
SLP(NELJSCH = SY(NEL))7DX 


a0 


cs SCAMeULdce Gne ON FECTIVE FUNCTION. . 00. ccccciccdivicccccscescce 
OBJ = 0.0 


DO 100 I=1,NEL 
BAVE(I) = (BASE(I)+BASE(I+1))72.0 
OBJ = OBJ + BAVE(I)XSQRT(1 + SLP(I)xx2) 
100 CONTINUE 
€ eee -CALCULATE ALPHA(I),BETA(I), AND BR EINGL ) 0 cicrs ERs aw ooo 
DO 150 I=1,NEL 
ALPHA(T JSATAN2(SLP(I )¥DX, DX) 
BETA(I)J=(PI72) = ALPHA(T) 
ELEN(I)J=SQRT(DX¥R2 + (SLP(I)XDX 2) 
150 CONTINUE 


---- determine the design Constraints. ....cccccscscccccecescece 
CALL ARCH_STRESS 


DO 230 I=1,NSNP 
G(I) = (SIGMA_LT(IJ/YIELD - 1.0) 
230 CONTINUE 
END 
C36 36 HE 3 3 3 IE 2 DE IE DE DE HE 3E DE HE DE 3 DE 3 CDE BE DE IE DE IE IE DE 9 IE DE IE IE 3K 3 IE 3 DE IE BE 3 BE IE 3 DE IE BE IE 3 3 9K DE 3E IE IE BE DE DE DE 9E IE IE IE BE DE 3E BE IE 
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FILE: ARCHOPT FORTRAN A 


SUBROUTINE ARCH_STRESS 


THIS SUBROUTINE IS USED TO PERFORM THE FINITE ELEMENT ANALYSIS 
OF THE STRESSES DEVELOPED IN AN ARCH OR BEAM FOR A GIVEN LOAD- 
ING. 


o«- Geclare the variables... ..cccccccccccscccccce cee tne ote 
INCLUDE "'ARCHCOM FORTRAN' 
INTEGER IPVT(99) 
REAL F(PG) 
REAL®8 BK(P4,P4),BF(P4),BU(P4),FAC(9801),WORK(99) 


qgg00q00n02 


-o--form the element and system matrices... ..c.ccecccccvccvvcces 
CALL FORM 


eo eee FORM THE FORCE VECTOR, PF icwieeeid cio. « caddiawiorecesokeususuebeeee aietete aes ete 
CALL FORCE_VECTOR (NEL,NDOF, ELEN, ALPHA,BETA, FA>F ) 


oes SET THE BOUNDARY CONDITIONS AND LOADS .« 00000 6 + Hees eee oe 
CALL BNDARY (NDOF,GK,CLAN,FX,FY,FM,F,BX1,BY1,BM1,BX2,BY2,BM2) 


eee SOLVE THE SYSTEM OF A So 56a 56 oo veo n0ko co cousiuaue ouciehedeusasue.e.s ene 
IF (PRCSN .EQ. 2) THEN 

-«+«eCHANGE GK AND F ARRAYS TO DOUBLE PRECISION........ 

CALL UPSCALE (NDOF,GK,F,PEK,BF) 

eee SOLVE THE SYSTEM OF EQUA TIGNS:. <icne.cccuste ouolt emaislcie aleve sieleterels steterenenene 
CALL DL2ARG (NDOF,BK,P4,BF,1,BU,FAC,IPVT, WORK ) 

e«-eCHANGE BU ARRAY TO SINGLE PRECISION... ccc ccccccccvcccscceecs 
CALL DOWNSCALE (NDOF,EU,U) 

EESE 

eee SOLVE THE SYSTEM OF EQUATIONS... ccc ccvcccrcccescceseccccccscccs 
Cee e ae (NDOF,GK,P4,F,3,U,FAC, IPVT, WORK ) 

ENDI 


qo0o06UCUCNTUCUCONOOCUCONCOTCUCONONLCUCONO 


oO 


CALL STRESS 


RETURN 
END 

C3 3 3 3 HE IEE 3 BE IE BE IE DE DE 3 BE DE DE IE DE BE BE DE DE HE IE DE DE DE BE 3 HE IE BE IE IE DE DE DE IE 3 DE DE DE IE DE DE DE DE DE DE DE 3 IE DE DE DE DE DE DE BR DE DE IE DE DE HE 3 
SUBROUTINE FORM 


Cc 
C «+e. determine the stress distribution. .....<««ccles ols ciecieicle aun 
Cc 
Cc 


Cc SH SB SSS SBS SSBB BS SBSS SSR EBB SS SSR BE SBS SRE SES SS BSB BSS S ES SS SaSS SSeS 2] 35S = 55 
c This subroutine is used to construct the global stiffness mat- 
c RIX FOR THE ARCH PROBLEM. 
C SS oe SS SSS SSS SS SBS SB SSB] S363 56344 63HSS6SSBSBBS6823B S82 38SSSS3e5 SSS S3$SSs== 
c -.- «Ge@Gtere thervariables......« « «<smeiesm- o © 6 yageie oceiiere lee <bete sence 
INCLUDE "ARCHCOM FORTRAN' 
INTEGER IEL,I,J,K,II,JJ,KK, III, JJ 
REAL C1,C2,C03,C04,C5,CA,CB, EK(P1,6,6)» GAMMA(6,6),EKGA(6,6)>, 
A : GAEKGA(6,6),BH»,BH3 
DO 120 IEL=1,NEL 
C ----Gefine the. comstarvts. ©xX. ccs dsc o- sieieieleisic susesebeiels + oe « cts 
Cl = YOUNG/ELEN(CIEL) 
C2 = (1.0/ELEN(IEL) )¥*%2.0 
C3 = (1.0)7(2. OXELEN(IEL)) 
C4 = (1.0)73.0 
C5 = (1.0)76.0 
Cc 
G eoeedimatialize the work arrays. « « oicsssis + ele o cheieeetetetsts) elcleieicia « « o seis 
DO 100 I = 1,6 
DO 90 J= 1,6 
EKPR(IEL,I,J)} = 0.0 
GAMMA(I,J) = 0.0 
EKGA(I,J) = 0.0 
GAEKGA(I,J} = 0.0 
EK(IEL,I,J) = 0.0 
90 CONTINUE 
100 CONTINUE 
c »..eCalculate the area and inertia termS....ccccccccccccccvece 
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FILE: ARCHOPT FORTRAN A 


BH = BAVE(TEL)XHGT 
BH3 = BAVE(TIEL )¥(HGT#3. 0) 


Cc 
Cc ~---- determine the EKPR PSRs 1 Motors coc ce cchaee es eetee eoeee#ee¢?e @eeeeseeeeve ee 
EKPR(IEL,1,1) = C1xBH 
EKPR(IEL,1,4) = -C1x%BH 
EKPR(IEL,2,2) = Ci1xC2xEBH3 
EKPR(IEL,2,3) = C1xC3xXBH3 
EKPR(IEL,2,5) = -C1xC2xBH3 
EKPR(IEL,2,6) = C1xC3*BH3 
EKPR(IEL,3,2) = C1xC3xXEBH3 
EKPR(IEL,3,3) = C1xXCGXBH3 
EKPR(IEL,3,5) = -C1xC3x¥BH3 
EKPR(IEL,3,6) = C1xC5SxXBH3 
EKPR(IEL,4,1) = -C1x*BH 
EKPR(IEL,4,4) = C1xEH 
EKPR(IEL,5,2) = -C1x¥C2xXBH3 
EKPR(IEL,5,3) = -C1xC3xBH3 
EKPR(IEL,5,5) = C1xC2xBH3 
EKPR(IEL,5,6) = -C1xC3xXBH3 
EKPR(IEL,6,2) = C1xC3xXEBH3 
EKPR(IEL,6,3) = C1xCSxXEHS3 
EKPR(IEL,6,5) = -C1xCSX¥BH3 
EKPR(IEL,6,6) = C1xCGxXBH3 
Cc 
(e: ~-- Getermine the GAMMA Segre tara eto eae Wot... ac laWaewalle le @ © 00 © €.0iexeve 6 
CA = COS(ALPHA(IEL) ) 
CS = COS(BETA(IEL)) 
GAMHA(1,1) = CA 
GAMMA(1,2) = CE 
GCAMMA(2,1) = -CB 
GAHHMA(2,2) = CA 
GAMHMA(3,3) = 1.0 
GAMMA(G4,4) = CA 
GAMMA(G,5) = CB 
GAMMA(5,4) = =-CB 
GAHMA(5,5) = CA 
GAMMA(6,6) = 1.0 
Cc 
c (emGmieeiMmitie CMe EKGA IARRAY. ..c cs cc ccc ccc cee cc cccc cc esec ese 
Po 220 I = 1,6 
BO 215 J = 1,6 
DO 210 K = 1,6 
EKGA(I,J) = EKGA(I,J) + EKPR(IEL,I,K)¥GAMMA(K,J) 
210 CONTINUE 
ZS CONTINUE 
' 220 CONTINUS 
Cc ee Getermine the GAEKGA APPAYV. ccc crcerecescececcesscccsccccce 


DO 240 I = 1,6 
DO 235 J = 1,6 
DO 230 K = 1,6 
GAEKGA(I,J) = GAEKGA(I,J)+GAMMA(K,I )¥EKGA(K, J) 
230 CONTINUE 


235 CONTINUE 
240 CONTINUE 
C 
c ~.--copy the GAEKGA array into the EK array... .ccccccccccccccecs 


DO 260 I = 1,6 
DO 250 J = 1,6 
EK(IEL,I,J) = GAEKGA(I,J) 


250 CONTINUE 
260 CONTINUE 
i 120 CONTINUE 
c iin ihitialize the GK APTAV cece ccecrcesecvresccesecresevesseece 


DO 150 I = 1, NDOF 
DO 140 J = 1, NDOF 
GK(I,J)} = 0.0 
140 CONTINUE 
150 CONTINUE 


Ue) 


FILE: ARCHOPT FORTRAN A 


Cc 
Cc eee Construct the GK Matrix. « ese 0c co aieelele ole ceric nennneen. ame 
DO 300 IEL = 1, NEL 
II = 3x¥(IEL-1) 
DO 290 J = 1, 6 
JJ = II + J 
DO 280 K = 1, 6 
= II + K 
GK(JJ,KK) = GK(JJ,KKIJtEKCIEL,J,K) 
280 CONTINUE 
290 CONTINUE 
300 CONTINUE 
Cc 
Cc RETURN 


END 
Ce Be BE HE HE HE HE HE HE HE HC HE HE IE WE EE HE OE BE BE UE IE 3 OE BE BE 3 BE HE BE OE BE EE I IE ME OE OE BE BE UE HE UE OE EE BE OE FE UE UE BE EM HE 3 OE OE IE OE OE EE 


SUBROUTINE FORCE_VECTOR (NEL,NDOF,ELEN,ALPHA,BETA,FA,>F) 


Cc This subroutine is used to construct the force vector for the 
c FEM PROBLEM SPECIFIED. 
Cc 


#0 ce DECLARE THE VARIABLES 6 o.c:e « aise 6 0 cierepabeneielsiels <jeketecs co 7cmeeeener nanan 
INTEGER NEL,NDOF,I,1I1,12,13,P1,P4% 


PARAMETER(P1=32,P4=99 ) 
REAL ELEN(P1),ALPHA(P1),BETA(P1),FA,F(P4) 


ORM THE PeNE COR. os eee ee ee 8 8 
= (ELEN(1)72.0) ® (-COS(BETA(1))) 
= (ELEN(1)72.0) % (CCSCALPHA(1))) 
= 0.0 


D9 100 I=2,NEL 
Il = (I-1)%3 + 1 
I2 = (i-1)*3 + 2 
= (i1-1)*3 + 3 


F(t) = (ELEN(I)72.0)"(-COS(BETA(I 
+(CELEN(I-1)72.0)*(-COS(EETAL 

I 

{ 


oeee@eeeteeeseeeeeeeeeeeeee @ 


))) 
1-1) 0 
y)) 
| 


“FCI2) = (ELENCI)72.0)*%(COS(ALPHA( 
+CELEN(I-1)72.0)%(COSCALPHA 
Ft = 0.0 
100 CONTINVE 


c 
F(NDOF-2) = (ELEN(NEL)/2.0)%(-COS(BETACNEL))) 
Tae = (ELEN(NEL)72.0)%(COSCALPHACNEL) )) 
F(NDOF = 


«eee SCALE THE F-VECTOR BY FAseccceae co ale cate cimiete s ciuka’stelebaneete’s clele 
DO 200 I=1,NDOF 

F(I) = FAXF(T) 
ac CONTINUE 


c RETURN 
END 
(C36 3 3 He Be 3 26 3 IE DE 3 BE DE 3 IE DE DE DE DE DE 3 IE DE DE 3 DE 3 ED IE 3 IE DE 3 DE DE 3 HE DE DE I 3 3 DE DE ED DE 3 3 IE DD 3 IE 3 3 DE 3 DE 3 9 3 2 3 3E 
pene ane eM ae »GK,CLAN, FX,FY,FM,F,BX1,BY1,BM1,BXe;, 
BY2,BM 


Cc 

[e: This subroutine is used to impose the boundary conditions upon 
Cc THE GLOBAL STIFFNESS MATRIX AND FORCE VECTOR. 
Cc 
Cc 


e 


a0 


~..-.-declare the VAP DSO LOS vcas 5. 0s 0: 0:00:00. 0:0 were: oo 4 lesa eeeeeeicus ia a tetenare 
INTEGER NDOF,BX1,BY1,BM1,BX2,BY2,BM2,CLAN,I,N,I1,12,13,P4 
PARAMETER(P4=99 ) 
REAL GK(P4,P4),FX,FY,FM,F(P4) 


c «eee invoke the essential boundary conditions.......cccccsccccce 
IF (BX1 .EQ. 1) THEN 
CALL IMPOSEBC (NDOF,GK,1,F) 
ENDIF 
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Cc 
IF (BY2 .EQ. 1) THEN 
CALL IMPOSEEC (NDOF,GK,2,F) 
ENDIF 
Cc 
IF (BM1 .EQ. 1) THEN 
CALL IMPOSEBC (NDOF,GK,3,F) 
ENDIF 
Cc 
IF (BX2 .EQ. 1) THEN 
N=NDOF-2 
CALL IMPOSEBC (NDOF,GK,N,F) 
ENDIF 
Cc 
IF (BY2 .EQ. 1) THEN 
N=NDOF-1 
CALL IMPOSEBC (NDOF,GK,N>,F) 
ENDIF 
c 


IF (BM2 .EQ. 1) THEN 
CALL INPOSEEC (NDOF,GK,NDOF,F) 
ENDIF 


c 

C ....-ADD THE CONCENTRATED LOAD TO THE FORCE VECTOR............ 
T1=(CLAN-1)x341 
T2=(CLAN-1 )x3+2 
IS=(CLAN=1 )¥34+3 


F(I1)J=FCI1)+FX 
F(I2)=F(I2)+FY 
FCI3)SFCI3)+FM 


Cc 

c RETURN 
END 

Ck De IE DE IE IE HE DE BE DE HE DE 3 3 DE IE 3 3 3 3 3 3 3 3 DE 3 3 3 EE 3 3 3 DE 3 3 3 3 3 3 3 De 3 0 3 3 3 3 De De 3 3 3 D3 3 3 3 3 3 IE IE IEE 
SUBROUTINE IMPOSEBC (NDOF,GK,N,F) 


G ss sea ssa caoaecoaa=casoaerecnaoaosaaoocoansaacoaonsaoaoansaoraaosoaoassens mamas noosa reas asaeaeneaana ae a 
c This subroutine is used to do the redundant leg work of impos- 
Cc ING THE BOUNDARY CONDITIONS. 

| Ot ee ee 
Cc eee e DECLARE THE VARIABLES.... eeeees? e@eeee?#e# eoeee@envnes8tcedsceeseeeee eee @ 


INTEGER NDOF,N,I,P4 
PARANETER(P4=99 ) 
REAL GK(P4,P4),F(P4) 


eee IMPOSE THE BOUNDARY CONDITION ON THE GK AND F ARRAYS...... 
DO 100 I=1,NDOF 
GK(N,I) = 0.0 
100 CONTINUE 
GK(N,N) = 1.0 


00 


F(N) = 0.0 
Cc 
Cc RETURN 

END 


C3 3 3 He 3 DE 3 DE 3 DE DE HE DE DE DE DE 3 DE DE DE 3 IE DE BE 3 DE DE IE BE DE IE DE 3 BE DE DE DE BE DE JE DE DE DE 3 DE DE 3 DE DE DE 3 DE 3 DE IE DE DE DE DE HE DE DE IE DE DE DE HE DE HE 3K 
SUBROUTINE UPSCALE(NDOF,GK,F,BK,BF) 


Cc 

c This subroutine is used to change the stiffness matrix & force 
Cc VECTOR FROM SINGLE PRECISION TO DOUBLE PRECISION IN ORDER TO 

Cc SOLVE THE LINEAR SYSTEM OF EQUATIONS IN DOUBLE PRECISION. 

Cc 
Cc 


eee +e DECLARE THE ESOS ES eS eee As C0 eer 
INTEGER NDOF,I,J,P4 

PARANETER (P4=99) 

REAL GK(P4,P4),F(P4) 

REAL®3 BK(P4,P4),BF(PG) 


c «+e GENERATE THE DOUBLE PRECISION COMPLIMENTS OF GK AND F...... 
DO 110 I=1,NDOF 
DO 100 J=1,NDOF 
BK(I,J) = GK(I,J) 
100 CONTINUE 


81 


FILE: ARCHOPT FORTRAN A 


BF(I) = F(T) 
110 CONTINUE 
Cc 
Cc RETURN 
END 
C3 BE IE HE 3 BE IE BE IE HE DE HE DE DE HE DEE DE EE J De ED EE 3 DD 3 3 3 3 3 UE EE DE DE HE IE IE 
SUBROUTINE DOWNSCALE(NDOF,BU,U) 


c 

c This subroutine is used to do down scale the double precision 

c SOLUTION OF THE LINEAR SYSTEM OF EQUATIONS BACK TO SINGLE PRE- 
c CISION. ADS COULD HAVE PROBLEMS WITH DOUBLE PRECISION NUMBERS! 
Cc 
€ 


-...declare the variables......ccccccccoee Te cee oe 
INTEGER NDOF,I,P4 
PARAMETER (P4=99) 
REAL UCPG) 
REALX8 BU(PG) 


c «.- GENERATE THE SINGLE PRECISION COMPLIMENT OF BU...... 
DO 100 I=1,NDOF 
Uli) = BUCi) 
100 CONTINUE 


Cc RETURN 
END 

Ck Ie HE I BE EI HEE HE DE BE HE DE DE HE HE DE DE DE HE I DE DE HE IE DE ED IE DEE DU DE DE HE DE UE DE DE HE DE IE DE DE DE DE DE DE DE DE DE HE DE DE HE DE DE HE DE HE HE HE 
SUEROUTINE STRESS 


Cc oe e eMOCLaratiOns. <<... cc cc ecco a8 0b ee 0c 6 6 ceeeneenere iene: 

INCLUDE 'ARCHCOM FORTRAN ' 

INTEGER T1,12,13,14,15,16,1,J5,K 

REAL CA1,CB1,K1,K2,FPR(P4,6),UPR(6),NORM1,NORM2, 
BEND1,BEND2 


-determine local forces from stiffness and displacement.... 
0 I=1, NEL 

T1=(I-1 )xS+1 

T2=(1-1)x3+2 

T3=(i-1)x%3+3 

I4=( i) 341 

I5=(1)*¥St2 

T6=(1)*¥3+3 


CS1= COS(BETA(T)) 
CA1= COS(CALPHA(I)) 


UCI1)*CA1 
-UCI1)*CB1 


wo) 


He ee 
By) 
7 
« 
Oo 
> 
= 
++ ++ 


~- FAFA 


250 CONTINUE 


DO 300 J=1,6 
DO 350 K=l, 
FPR(I,J) 
350 CONTINUE 
300 CONTINUE 
100 CONTINUE 
Cc -...-determine the bending and normal stresses... cccccccccccccs 
SIGMA_N(1) = ABSCFPR(1,1)%(1.0/(BASEC1 )¥HGT))) 
SIGMA _B(1) = ABSCFPR(1,3)%(6.0/(BASE(1 )¥(HGTX®2.0)))) 


6 
= FPR(I,J) + EKPR(I,J,K)¥UPR(K) 


SIGMA_T(1) = SIGMA_B(1) + SIGMA_N(1) 
DO 400 i=2,NEL 

Kl = 1.0/(BASECI )*HGT) 

K2 = 6.0/(BASE(I )¥(HGTx#2.0)) 
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NORM1 = ABS(FPR(1,1)%K1) 

NORM2 = ABS(FPR(1i-1,4)%K1) 

BEND1 = ABS(FPR(i,3)®K2) 

BEND2 = ABS(FPR(i-1,6)*K2) 
SIGMA_N(i) = (NORM1+NORM2)72.0 
SIGMA_B(iJ) = (BENDI+B5END2)72.0 
SIGMA_T(i) = SIGMA_B(i) + SIGMA_N(i) 


400 CONTINUE 
SIGMA_N(NSNP) = ABS(FPR(NEL,4)¥(1.0/(BASE(NSNP )¥HGT ) )) 


SIGMA_B(NSNP } ABS(FPR(NEL > 6) 
: (6.0/7 (BASE(NSNP J¥(HGT¥%2.0)))) 
SIGMA_T(NSNP) = SIGMA_B(NSNP) + SIGMA_N(NSNP) 
c 
c RETURN 


END 
3K 3K HE BE BE HE HC DE I UE BE SE BE DE BE UE UE UE IE HE IE I TE 3 HE 3 DE 3 BE SE UE DE SE SE SE DE SE IE SE IE FE IE IE OE FE IE SE DE IE OE UE OE DE OE DE FE DE DE DE OE OE 3G DE OE DE DE OE OE OE 


SUBROUTINE ARCH_OUTPUT 


c 

Cc THIS SUBROUTINE FORMATS THE FINAL RESULTS AND OUTPUT OF THE 
Cc OPTIMIZATION PROBLEM AND STORES IT IN A FILE NANED ARCH_OUT 
C 
c 


eee eae = SBS @ew ees eweese ese ee @ ee @ SOBs Os FOO ZB SF fF SB SB SF HB SO BPS es OF SOF SF ZT SF ZT es TT ee ee ese es Se |S se 
Se @@@e2 Ss Ss Oe SB SS eB BS SS SF SF SS SF SBS SF SF SSBF SF SF SBF SF SSF SF SF SF SSF FZ Ss Ze TF Fs 8s ses SF SS SS SO SF SF eee eS Se 


ee.» DECLARE V/AIRFCA En tSterens oxoue oo slenenelelshele ¢ 6 ekene 66 ee ete eens. « eeees? e@eee 
INCLUDE "ARCHCOM FORTRAN‘ 
REAL VOL, VOLUME 


«eee OPEN OUTPUT FILE AND WRITE POPE ss oleleieie « 6 60eleie0e« «< eaeoe2ees? & 
OPEN(9, FILE='ARCHOUT', STATUS='OLD' ) 


WRITE(9,100) LABEL 

WRITE(9,100) ° OPTIMIZATION SOLUTION’ 

WRITE(9, 105) ' wqwe mm n mw nn enn wenn ene ewe wwe ene nen nee n nen ne ene ne 
100 FORMAT(75X,A} 
105 FORMAT(SX,A} 


c 
Cc veo SEC TPON ER een UC. cag gig alias. « o's 0 68S 6 4.018 siereneee.e e@eee#e e 
WRITE(9,100) ' AJ PROBLEM PARAMETERS: ' 
WRITE(9,110) * HORIZ. SPAN:', L, ' YOUNGS MODULUS: ', YOUNG 
WRITE(9,110) ' VERT. SPAN :', H, ' YIELD STRENGTH: ', YIELD 
WRITE(9,115) ' NO OF DESIGN VAR:', NDV, ' NO OF ELEMENTS:',NEL 
110 FORHAT(8X,A,F12.3,T38,A,F12.1) 
115 FORHAT(8X,A,17,T38,A,110) 


c 
Cc ieee SERGI LON gh I eoeeeee @eeteoeeeeeeeoeeeeeeeeeeeeee ee eee eG @ 
WRITE(9,100) ' B) DERIVED CONSTANTS: ! 
WRITE(9,120) ' NO OF SYSTEM NODAL POINTS...‘ »NSNP 
WRITE(9,120) * NO OF DEGREES OF FREEDOM....',NDOF 
WRITE(9,125) * HORIZ. LENGTH PER ELEMENT...',DX 
120 FORMAT(8X,A,16) 
125 FORMAT(8X,A,F12.4) 


a0 


Cc 

Cc eee -Ssection SEO IO” AI ee 
WRITE(9,100) ' ©) STRUCTURE LOADING: ' 
WRITE(9,125) EREOTEET TS ccs oc slele oe eters acco ecee FFX 
WRITE(9,125) SINE ET oo) soe reels cleo os 006 oe PY 
WRITE(9,125) SETAM PMNS Gs 5 7stp ko «as ass ss ies oe SEM 
WRITE(9,125) SN cna eeceoee cesses esesesceee OFA 
WRITE(9,120) ‘CONCENTRATED LOAD AT NODE...',CLAN 


eee -SeCtion WA pc ooo oa WW 6 A 6 6 eee 
$X(1)=0.0 
S¥(1)=0.0 
DO 150 I=1,NEL 
SX(I+1)=Ix¥DX 
SYCI+1)=DX¥SLP(I) + SY(T) 
150 CONTINUE 


ok») 


é3 


wv 


FILE: 


a0 


340 


400 
410 
420 
430 
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WRITE(9,100) * D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: ' 
WRITE(9,210) ‘NODE’, 'X-COORD', "Y-COORD','BASE', "AREA' 


FORMAT(8X,A,T21,A,T36,A,T49,A,T62,A) 
FORMAT eke tis TL T+ FlGasal 32> F 10-5) TS Sa Raia aaaniaal 
VOLUME = 0. 


DO 300 I=1,NSNP 

AREA = BASE(I)*HGT 

WRITE(9,220) I,SX(I),SY(I),BASE(I),AREA 
CONTINUE 


~--- section MEU, occ c 6 oes 6 Dae eee e@o@eee8 @ eoeeeeeeee eee ee ee ee @ @ 
VOL = HGT¥DXXOBJ 
WRITE(9,100) * E) OBJECTIVE FUNCTION: ' 


WRITE(9,310) " TOTAL STRUCTURE VOLUME: ',VOL 
FORMAT(712X,A,F12.67) 


WRITE(9,330) "NODE','NORMAL STRESS', "BENDING STRESS','TOTAL' 
DO 320 I=1,NSNP 

WRITE(9,340) i,SIGMA_N(i),SIGMA_B(i),SIGMA_T(i) 
CONTINUE 
FORMAT(8X,A,T18,A,T35,A,T57yA) 
FORMAT(8X,14,T15,F14.1,132,F14.1,T49,F14.1) 


sine SCCCLOND "Rien. sos aus eee Oe on Sereno tS 
WRITE(9,100) ' F) BOUNDARY CONDITIONS: ' 
WRITE(9,410) "NODE','X-DISPL','Y-DISPL', SLOPE’ 
WRITE(9,430) 1,EX1,BY1,BM1 
WRITE(9,430) NEL41,BX2,BY2,BM2 


ooo eo SECTION MGM, 2 oo 0.0 2/5 o otuMW slctmeuec.sicisle © 0 alsseceisieleisipisisiuiekessl> + « « ocolele 
WRITE(9,100) ° G) SOLUTION VECTOR: ! 
WRITE(9,410) ‘NODE','X-DISPL','Y-DISPL', "SLOPE' 
DO 400 I=1,NSNP 

T1=(I-1)x%3+1 

T2=(I-1)x*3+2 

I3=(I-1)x*3+3 
WRITE(9,420) I,U(T1),U(T2),U(T3) 
CONTINUE 
FORMAT (T9,A>T17,A,T31,A,T455A) 
FORMAT(7X,15,3E14.6) 
FORMAT(7X,1I5,T20,14,T34,14,T48,14) 


RETURN 
ED 
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OPTIMIZATION NO.1 
OPTIMIZATION SOLUTION 


A} PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: le 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDCM.... “Ie 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
Bees. ss eeoee7eeee#eee @eeeoeeteese7ee 0.0000 
PY ears en teteva@ 6), 6's oe eles eee e @ 406 @ 0.0000 
BES a a oie ete cae 6: ace Stee erelec ale 6 6 © © 0.0000 
PPA iaic iat ckeussleuanecacateyetaVelewayatesersie 0 « -100.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 0.15000 0.22500 
2 2.66667 4.30175 0.15000 0.22500 
S S2esss3 8.23110 0.15084 Ovecca/ 
& 8.00000 11.74150 0.15022 Oe2e533 
> 10.66666 14.87895 0.15092 0.22638 
6 l3eo55>5 17.72124 0.15014 Gmecsce 
7 16.00000 20-31323 0.15000 0.225090 
8 18.66666 22.73409 0.15000 0.22500 
9 21.35533 24.91997 0.15000 0.22500 
10 23.99998 26. 86490 0.15000 0.22500 
11 26.66666 28.66069 0.15000 0.22500 
le 29.33331 30.59526 0.15000 0.22500 
13 31.99998 31.99998 0.15000 0.22500 
E) OBJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 10.324365 
NODE NORMAL STRESS BENDING STRESS TOTAL 
1 3$6039.3 O.1 36039.5 
e 36031.7 T5159.3 51191.0 
3 goses.3 pol75.0 51998.3 
& 35775.3 P5313.8 51289.2 
2) 35809.8 15699.2 a1509.1 
6 35994.7 16009.5 52004.3 
7 36033.5 15762.4 51796.0 
8 36044.7 9766.7 45811.4% 
9 36049.8 7516.0 43565.8 
10 36041.4 12457.1 48498.6 
it 36016.8 16026.7 52043.5 
12 36039.7 15969.3 52009.0 
13 36081.4% 0.0 36081.4% 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
1 1 Z 0 
13 i 1 0 
G) SOLUTION VECTOR: 
NO X-DISPL Y-DISPL SLOPE 


Yo 
NOU D Wem 


0.509330E-01 
0.863026E-01 
0.106222E+00 
0.114465E+00 
0.113685E+00 
0.105703E+00 


0.000000E+00 0.000000E+00 
-0.387271E-01 
-0.696029E-01 
-0.913596E-01 
-0.104831E+00 
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-0.131531E-01 
-0.114481E-01 
-0.814113E-02 
-0.505612E-02 
-0.221819E-02 
-0.110495E+00 0.528031E-03 
-0.108688E+00 0.315360E-02 


CASE1 A 


0.921987E-01 
0.761242E-01 
0.589296E-01 
0.396341E-01 
0.162884E-01 
0. 000000E+00 


~0.100249E+00 
~0.871751E-01 
-0.703301E-01 
-0.485910E-01 
-0.231430E-01 

0.000000E+00 


0.519679E-02 
0.652113E-02 
0.798611E-02 
0.100210E-01 
0.100252E-01 
0.895565E-02 
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OPTIMIZATION NO.1A 
OPTIMIZATION SOLUTION 


A) PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: 2 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 3 
NO OF DEGREES OF FREEDOM.... =o 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
acter lccls c/s s ses 6 60 wile sles o 0 60 6 0 0.0000 
FY spereweeeneer ss: . . SUE. ci ccc cee 0.0000 
Tels © ones ss sues esisevecc sv ces 0.0000 
2 Cee ee -100.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 0.79198 1.18797 
2 2.66667 2.66933 0.78408 1.17612 
3 5.33333 5.33867 0.78335 1.17503 
& 8.00000 8.00800 0.98644 1.47966 
5 10.66666 10.67200 1.16866 1.75299 
6 13.33333 13.33852 1.27656 1.91424 
7 16.00000 16.00357 1.30861 1.96291 
8 18.66666 18.66756 1.27457 1.91186 
=) 21.33333 21.33154 1.16462 1.74693 
10 23.99998 23.99666 0.98731 1.48096 
1d 26.66666 26.66280 0.78324 1.17486 
IZ 29. s5seu 29.32883 0.78168 1.i7252 
13 31.99998 31.99998 0.78982 1.18473 
E) OSJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 68.386749 
NODE NORMAL STRESS BENDING STRESS TOTAL 
1 91.5 0.0 oi. 5 
2 92.6 26676.4% 26769.0 
3 92.2 48556.6 48648.8 
4 1a.9 52068.1 52140.9 
s 61.6 52094.5 S2iso. 1 
6 56.0 52ie3.3 52239.4% 
r 54.3 52353.9 52408.3 
8 56.5 52263.8 52320.3 
9 62.2 52273.6 52eco-8 
10 73.0 52014.1 52087.1 
Zl 92.3 48564.5 48656.7 
12 91.5 26770.9 26862.4% 
13 89.5 0.2 89.8 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
1 1 i 0 
13 a 1 0 
G) SOLUTION VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 
1 0.000000E+00 0.000000E+00 -0.432361E-01 
2 0.113423E+00 -0.113326E+00 -0.410106E-01 
3 0.215283E+00 -0.215100E+00 -0.347033E-01 
G 0.297166E+00 -0.296917E+00 -0.262325E-01 
5 0.355750E+00 -0.355571E+00 -0.17507GE-01 
6 0.390952E+00 -0.390785E+00 ~-0.876822E~-02 
7 0.402699E+00 -0.402549E+00 -0.976327E-05 
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FILE: 


FILE 


CASELIA A 


0.390992E+00 
0.355829E+00 
0.297182E+00 
0.215351E+00 
0.113571E+00 
0. 000000E+00 


-0.390840E+00 
-0.355653E+00 
-0.296983E+00 
-0.215151E+00 
-0.113364E+00 
0. 000000E+00 


0.875361E-02 
0.175100E-01 
0.262491E-01 
0.347106E-01 
0.410251E-01 
0.432590E-01 
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FILE: FICE CASE2 A 


OPTIMIZATION NO.2 
OPTIMIZATION SOLUTION 


A} PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: l2 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C} STRUCTURE LOADING: 
F eeeee eoeeoeeee ee OO ee Gee ee Oe Oo 0.0000 
0 RRP Saar ar sey ace -2000.0000 
Faliliclco © oo 6 0s e oo 90s si vicecccecse 0.0000 
Fpeieielelsis es 2) = © soe s+ <iei © seis 0.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 0.15000 0.22500 
2 2.66667 3.00301 0.15000 0.22500 
3 5.33333 6.30157 0.15000 0.22500 
G 8.00000 9.69067 0.15000 0.22500 
2 10.66665 12.91733 0.15000 0.22500 
6 13. 33333 16.27870 0.15000 0.22500 
U 16.00000 19.31017 0.15000 0.22500 
8 18.66666 21.33359 0.15000 0.22500 
7 21.33333 23.39781 0.15000 0.22500 
10 23.99998 25.46082 0.15000 0.22500 
ual 26. 66666 27.54510 0.15000 0.22500 
l2 29.33331 29. 74944 0.15000 0.22500 
13 31.99998 31.99998 0.15000 0.22500 
E) OBJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 10. 239457 
NODE NORMAL STRESS BENDING STRESS TOTAL 
i 31082.2 0.3 31082.6 
2 31095.7 20958.1 52053.9 
3 31109.7 18671.2 49780.9 
4 31109.4 16336.5 47445.9 
5 31107.6 12623.3 43730.9 
6 31097.4 5386.4 36483.7 
7 27975.7 24094.0 52069.8 
8 24863.7 27275.3 52139.0 
9 24863.7 27252.7 52116.4 
10 24863.8 27320.7 52184.5 
11 24857.5 25720.0 5O0S77.5 
le 24846.4 14678.1 39524.5 
13 24841.8 0.0 24841.8 
&) weoe CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
1 1 0 
13 1 1 0 
G) NGoE’ ne VECTOR: 
ODE X-DISPL Y-DISPL SLOPE 
1 0.000000E+00 0.000000E+00 -0.169057E-01 
2 0.461328E-01 -0.465305E-01 -0.150352E-01 
3 0.866829E-01 -0.849686E-01 -0.112998E-01 
G Q.115632E+00 -0.114025E+00 -0.799962E-02 
5 0.134663E+00 -0.134988E+00 -0.526077E-02 
6 0.146308E+00 -0.149906E+00 -0.354366E-02 
7 Q.151127E+00 -0.159717E+00 -0.898687E-03 
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FILE: 


FILE 


CASE2 A 


0.146949E+00 
0.134488E+00 
0.113602E+00 
0.840616E-01 
0.448603E-01 
0.000000E+00 


-0.158800E+00 
-~0.147268E+00 
~0.124838E+00 
-0.915992E-01 
-0.486744E-01 

0.000000E+00 


0.292259E-02 
0.700891E-02 
0.110977E-01 
0.150871E-01 
0.181930E-01 
0.193311E-01 
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Pree: PICe CASES A 
OPTIMIZATION NO.3 
OPTIMIZATION SOLUTION 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: 12 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
oS oe a Re ae ee 2000.0000 
BERING, coca aac. a> a: Wabelelateterelacs 6 4 se5e-* ~2000.0000 
Merits so 2c nisictevere © 6S eieicie's 6 eels 0.0000 
Pei e i ra S66) EEE es « 0.0000 
CONCENTRATED LOAD AT NODE... 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 0.15030 0.22545 
ec 2.66867 3.57692 0.15037 0122556 
3 5.33333 7.11881 0.15045 OS2e568 
4 8.00000 11.13781 0.15079 0.22618 
5 10.66666 14.64122 0.15017 0.22526 
6 13.33333 18.30432 0.15000 0.22500 
7 16.00000 21.47563 0.23559 0.35339 
8 18.66666 23.44431 0.15046 0.22568 
9 2loocess 25.24850 0.15008 0.22512 
10 23.99998 26.98657 0.15068 0.22603 
a1 26.66666 28.49719 0.15040 0.22561 
le 29.33331 30.15172 0.15009 0.22514 
13 31.99998 31.99998 0.15012 0.22517 
E) OSJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 10.841561 
NODE NORMAL STRESS BENDING STRESS TOTAL 
1 32852.1 0. 32852.2 
2 32835.8 8291.1 41126.9 
3 32805.9 19268.1 52074.0 
4 32730.9 6463.6 39194.5 
5 32878.3 7478.3 40356.5 
6 32873. 2 9143.1 42021.4% 
7 21034.9 31050.3 52085.2 
8 S33hi5.7 16387.2 49562.9 
9 33301.1 2650.6 35951.7 
10 33157.7 3055.2 36812.9 
11 33220.8 15609.0 48829.8 
ie 33293.2 18737.4 52030.6 
13 33272.5 . 33272.6 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
1 1 1 0 
13 1 1 0 
G) SOLUTION VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 
1 0.Q000000E+00 0.000000E+00 -0.927029E-02 
2 0.292592E-01 -0.279060E-01 -0.844808E-02 
3 0.520934E-01 -0.511700E-01 -0.573258E-02 
G 0.676226E-01 -0.677945E-01 -0.436322E-02 
5 0.806200E-01 -0.837384E-01 -0.426672E-02 
6 0.903615E-01 -0.969766E-01 -0.259323E-02 
7 0.912367E-01 -0.102322E+00 0.155546E-02 
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FILE: 


FILE 


CASE3 A 


0.815972E-01 
0.675751E-01 
0.529933E-01 
0.398219E-01 
0.227118E-01 
0. 000000E+00 


-0.940706E-01 
-0.797142E-01 
-0.638010E-01 
-0.474255E-01 
-0.264506E-01 
0. 000000E+00 


0.528782E-02 
0.665123E-02 
0.657931E-02 
0.739218E-02 
0.978718E-02 
0.111381E-01 
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FILE: FILE CASES A 


OPTIMIZATION NO.4 
OPTIMIZATION SOLUTION 


A) PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: 12 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
SOO DOO OOOHOS UAC DOO Cre 0.0000 
FN isn sess sc ces ee secsssccscecs 0.0000 
Clio 6 6 oe ere ea 6000.0000 
FA. oe@eeee ee @ @ 0.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
S) 0.00000 0.00000 0.15087 0.22631 
2 2.66667 2.69333 0.15000 0.22500 
3 5.33333 5.38662 0.15000 0.22500 
| 8.00000 8.07996 0.15000 0.22500 
iB, 10.666€6 10.77329 0.15000 0.22500 
6 13.33333 13.46662 0.15000 0.22500 
7 16.00000 16.10661 0.15377 0.23066 
8 18.66666 18.74661 0.15000 0.22500 
9 21.33333 21.43994 0.15014 0.22521 
10 23.99998 24.13327 0.15005 0.22507 
11 25.66666 26.77327 0.15000 0.22500 
12 29. 33331 29.41327 0.15000 0.22500 
13 31.99998 31.99998 0.15000 0.22500 
E) OEJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 10.207288 
NODE NORMAL STRESS BENDING STRESS TOTAL 
1 18.9 0.0 18.9 
2 19.0 8935.8 8954.8 
3 19.0 17871.6 17890.6 
G 19.0 26806.6 26825.6 
= 19.1 35742.1 35761.2 
6 22.0 44677.8 44699.8 
7 24.3 52024.6 52048.9 
8 22.0 44303.3 44325.2 
7 19.0 35335.9 35354.9 
10 21.8 26425.4% 26447.2 
it 24.6 l7577. I 17615.7 
12 27.6 8749.7 8777.3 
13 30.6 0.1 30.7 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
i 1 1 0 
13 1 1 0 
G) SOLUTION VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 


0.000000E+00 0.000000E+00 -0.911602E-02 
0.23880GE-01 -0.236406E-01 -0.836558E-02 
0.43710GE-01 -0.432711E-01 -0.610774E-02 
0.554325E-01 -0.548738E-01 -0.234471E-02 
0.549926E-01 -0.544349E-01 0.292347E-02 
0.383368E-01 -0.379406E-01 0.969686E-02 
0.238511E-02 -0.162140E-02 0.177837E-01 


SIO UID GOIN) 


hes 


FILE: 


FLEE 


8 
9 
10 
11 
12 
13 


CASES A 


-0.336480E-01 
-0.505593E-01 
-0.514326E-01 
-0.405827E-01 
-0.219884E-01 

0.000000E+00 


0.347801E-01 
0.515273E-01 
0.523953E-01 
0.414402E-01 
0.226625E-01 
0. 000000E+00 


0.975857E-02 
0.305128E-02 
-0.215084E-02 
-0.582135E-02 
-0.801783E-02 
-0.874020E-02 
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FILE: FILE CASES A 
OPTIMIZATION NO.5 
OPTIMIZATION SOLUTION 


A) PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: 12 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
Pieter erie also oie velels «cUMe oco.s.0 6 6 6 0.0000 
FevEV ea 506 6 4.6 6 's'.6 6.0 6 6 SeSiee o 6 6 0s 0 0.0000 
lle sce cc ese cescsescccsscses 0.0000 
Pee aici ian aekel oars lake Oks one dauelee © * -100.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 0.15000 0.22500 
2 2.66667 10.66666 0.16746 0.25118 
3 5.33333 16.43103 0.15000 0.22500 
4 8.00000 20.18451 0.15000 0.22500 
S) 10.66666 22.93553 0.15000 0.22500 
6 13. 33333 25.15350 0.15006 0.22510 
7 16.00000 26.92931 0.15000 0.22500 
8 18.66666 28.34477 0.15000 0.22500 
2 21.33333 29.47061 0.15090 0.22500 
10 23.99998 30.29001 0.15000 0.22500 
ll 26.66666 30.90901 0.15000 0.22500 
l2 29.33331 31.43495 0.15000 0.22500 
13 31.99998 31.99998 0.15000 0.22500 
E) OBJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 11.328694 
NODE NORMAL STRESS BENDING STRESS TOTAL 
B: 13798.3 0.1 13798.5 
2 12371.2 39641.7 52012.9 
3 13857.5 38158.4 52015.9 
4 13897.5 36095.7 49993.2 
5 13901.5 38115.7 52017.2 
6 13895.2 38108.5 52003.6 
7 13903.2 38094.6 51997.8 
8 13906.0 38087.9 51993.8 
9 13909.4 37119.4 51028.8 
10 13912.9 37831.8 51744.7 
11 13924.4 35752.2 49676.6 
12 13922.3 2o0ep.1 39477.4 
13 13910.7 1.5 13912.2 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
1 0 1 0 
13 1 1 0 
G) SOLUTION VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 
1 -0.595168E+00 0.000000E+00 -0.429371E-01 
2 -0.174663E+00 -0.110052E+00 -0.327189E-01 
3 -0.196828E-01 -0.184795E+00 -0.217266E-01 
G 0.462421E-01 -0.234246E+00 -0.141289E-01 
5 0.752630E-01 -0.264850E+00 -0.781040E-02 
6 0.848358E-01 -0.278872E+00 -0.193522E-02 
7 0.822190E-01 -0.277620E+00 0.349024E-02 


2 


FILE: 


FILE 


CASES A 


0.724252E-01 
0.587701E-01 
0.4462G41E-01 
0.307934E-01 
0.166607E-01 
0.000000E+00 


-0.262154E+00 
~0.233261E+00 
-0.191628E+00 
-0.137662E+00 
-0.725301E-01 

0.000000E+00 


0.860142E-02 
0.134390E-01 
0.180855E-01 
0.225620E-01 
0.262650E-01 
0.278131E-01 
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FILE: FILE CASE6 A 


OPTIMIZATION NO.6 
OPTIMIZATION SOLUTION 


A) PROBLEM PARAMETERS: 


NSOUDGN em 


-0.419812E-01 
0.630188E-01 
0.662482E-01 
0.522159E-01 
0.444548E-01 
0.407255E-01 


-0.792292E-01 
-0.118314E+00 
-0.122027E+00 
-0.112946E+00 
-0.104995E+00 
-0.990568E=-01 


-0.223869E-01 
-0.631173E-02 
0.299570E-02 
0.399436E-02 
0.258829E-02 
0.268114E-02 


o7 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN ; 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: l2 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
Bera cree arcs 0s 6.6. s 6.060%. 0 CEM 6 eo sce 0.0000 
aVeeeevetctare clave oats 6 ecsrele.e 6 eevetenacs 0.0000 
FM ierotevs vers ceise eeee3#2ee#8e?se# eoeoeoetee 8 @ @ 0.0000 
REET eteic ils 6 se .ls.cik 6s sieletslele ss -100.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
HODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 0.15000 0.22500 
2 2.66667 6.78236 0.30329 0.45494 
3 5.35555 14.31105 0.24584 0.36876 
4 8.00000 19.67329 0.15000 0.22500 
5 10.66666 22.94662 0.15000 0.22500 
6 13. So0o5 24.93712 0.15000 0.22500 
q 16.00000 26.01016 0.15000 0.22500 
8 18.66666 27.14098 0.15000 0.22500 
9 21.33sse 28.00891 0.15000 0.22500 
10 23.99998 28.76320 0.15000 0.22500 
11 26.66666 29.50047 0.15000 0.22500 
r2 29.33331 30.30164 0.15004 0.22506 
13 31.99998 31.99998 0.43413 0.65119 
E) OSJECTIVE FUNCTION: 
% TOTAL STRUCTURE VOLUME: 14.430920 
NODE NORMAL STRESS PENDING STRESS TOTAL 
1 12224.7 O71 12224.8 
2 6057.4 45945.8 52003.1 
3 7547.9 44G69.5 52017.4% 
G 12709.8 19397.3 32107.1 
5 13071.6 8753.1 21824.7 
6 13083.4 10261.1 23344.6 
7 12983.4 11714.5 24698.0 
8 12980.5 16215.1 29195.6 
9 12969.5 18953.1 31922.6 
10 12975.1 13532.8 26507.9 
11 12951.8 G652.1 17603.9 
2 12164.7 39871.4 52036.1 
13 3942.6 48101.8 52044.4 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
1 0 1 0 
is 1 1 1 
G) SOLUTION VECTOR: 
NOD X-DISPL Y-DISPL SLOPE 
-0.238121E+00 0.000000E+00 -0.323442E-01 


— tt 
ANH O0 0 


CASE6 A 


0.355767E-01 
0.296101E-01 
0.226288E-01 
0.147434E-01 
0.702276E-02 
0. 000000E+00 


-0.901286E-01 
-0.757135E-01 
-0.554353E-01 
-0.314062E-01 
-0.987708E-02 

0.000000E+00 


0.447894E-02 
0.667067E-02 
0.867125E-02 
0.921718E-02 
0.646195E-02 
0.000000E+00 
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FILE: FILE CASE7 A 


OPTIMIZATION NO.7 
OPTIMIZATION SOLUTION 


ees @ SB SB SO SB SB SO SB BS Se FZ 2 SSG es ee BSS FB SSB SF SB SB SF SF SF SB SF SBS STF FF SF TB SST SBS SF SF SB SO SF SO SB SB SS SFwZeEe SO SOF SF BF SB SF ee | ws 


A) PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: ié 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
oe ee a es 0.0000 
ON oo rr a. een 0.0000 
Eo 0 OS allel ae. ca ae 0.0000 
BEIGE! «06 wisup ss oo see ccs as cele -~100.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 1.31120 3.93359 
2 2.66667 2.66377 1.10192 3.30575 
3 5.33333 5.25856 0.91330 2.73989 
G 8.00000 7.97715 0.73940 2.21821 
5 10.66666 11.05349 0.57427 1.72280 
6 13. 33333 14.02065 0.43146 1.29439 
7 16.00000 16.77490 0.31308 0.93925 
8 18.66666 19.50882 0.30000 0.90000 
9 21. 33333 22.17821 0.30000 0.90000 
10 23.99998 24.47510 0.30000 0.90000 
11 26.66666 27.00243 0.30000 0.90000 
le 29. 33331 29.40999 0.30000 0.90000 
13 31.99998 31.99998 0.30000 0.90000 
E) O8JECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 72.614471 
NODE NORMAL STRESS BENDING STRESS TOTAL 
1 0.5 52070.8 52071.3 
2 S29 52066.6 52075.5 
3 16.1 52056.2 5207252 
G 61.1 51982.5 52043.6 
5 123.9 51879.7 52003.6 
6 107.4 51886.1 51993.4% 
7 92.3 51919.3 52011.6 
8 79.6 37074.2 37153.8 
9 60.3 23343.8 23404.1 
10 28.9 13394.2 13423.0 
ll 9.9 5932.8 5942.7 
12 7.7 1530.9 1538.6 
13 0.0 5.7 5. 
F} BOUNDARY ee. 
NODE X-DISPL Y-DISPL SLOPE 
1 1 1 1 
13 0 0 0 
G) SOLUTION VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.597678E-02 -0.598316E-02 -0.436127E-02 
3 0.230543E-01 -0.23530GE-01 -0.866592E-02 
G 0.528078E-01 -0.527179E-01 -0.13068GE-01 
5 0.100533E+00 -0.9G1091E-01 -0.177673E-01 
6 0.160387E+00 -0.147924E+00 -0.223668E-01 
7 0.228394E+00 -0.213783E+00 -0.267883E-01 
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FILE: FILE 


CASE7 A 


0.307128E+00 
0.392384E+00 
0.470241E+00 
0.558680E+00 
0.644186E+00 
0.736590E+00 


-0.290595E+00 
-0.375776E+00 
-0.466157E+00 
-0.559472E+00 
-0.654177E+00 
-0.749316E+00 


-0.305781E-01 
-0.331111E-01 
-0.345477E-01 
-0.353367E-01 
-0.356349E-01 
-0.356983E-01 
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FILE: FILE CASE8 A 


OPTIMIZATION NO.8 
OPTIMIZATION SOLUTION 


A) PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN : 32.000 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: 12 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
PiteCMMCM GPa cies ses oc oc c's 6 6500 5 6 6 « 0.0000 
Pi treet cia) alaleie'e) clots ehallencrelevece 6 eels ane 0.0000 
iMieke sos 606.6 ocolelere 6 600 @ 6 6 «se 86 0.0000 
A sas 8 ieee a as cose aialoceisceuauaneueres « -100.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD EASE AREA 
1 0.00000 0.00000 0.40294 0.60440 
2 2.66667 3.27370 0.15000 0.22500 
3 5.33333 8.18358 0.15000 0.22500 
4 8.00000 12.49695 0.15000 0.22500 
5 10.66666 15.68683 0.15000 0.22500 
6 13.33333 18.42828 0.15000 0.22500 
7 16.00000 20.90140 0.15000 0.22500 
8 18.66666 23.20117 0.15000 0.22500 
9 21.33333 25.28195 0.15000 0.22500 
10 23.99998 27.28166 0.15000 0.22500 
11 26.66666 29.15797 0.15000 0.22500 
12 29.33331 30.54808 0.15000 0.22500 
13 31.99998 31.99998 0.21518 0.32277 
E) OBJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 11.308629 
NODE NORMAL STRESS BENDING STRESS TOTAL 
st 9208.7 42801.0 52009.6 
2 25084.6 3375.4% 28460.0 
3 25443.1 6676.3 32119.3 
4 25460.7 9375.9 34836.6 
5 25436.8 3831.6 29268.3 
6 25386.7 18147.5 43534.2 
7 25363.4 26649.6 52012.9 
8 25365.1 25377.8 50742.9 
9 25379.9 20239.0 45618.8 
10 25414.8 1767.0 27181.8 
11 25514.8 25373.4% 50888.2 
T2 25579.6 26397.0 51976.7 
13 17825.3 34401.1 52226.4% 
F} BOUNDARY CONDITIONS: 
NODE X=-DISPL Y-DISPL SLOPE 
1 1 1 1 
13 1 1 1 
G) wae eg VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 
1 0.000000E+00 O.QQ00000E+00 0.000000E+00 
2 Q0.113937E=-01 -0.117173E-01 -0.568126E-02 
3 0.342988E=-01 -0.295478E=-01 -0.443323E-02 
G 0.505141E-01 -0.4466311E-01 -0.47374GE-02 
5 0.648279E-01 -0.611976E-01 -0.524971E-02 
6 0.749570E-01 -0.755686E-01 -0.338176E-02 
7 O.768717E-01 -0.821558E-01 0.238845E-03 
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0.693485E-01 
0.544237E-01 
0.346055E-01 
0.157549E-01 
0.522320E-02 
0. 000000E+00 


-0.779902E-01 
-0.635129E-01 
-0.417868E-01 
-0.198010E-01 
-0.514698E-02 

0.000000E+00 


0.431015E-02 
0.773890E-02 
0.936886E-02 
0. 765835E-02 
0.419866E-02 
0.000000E+00 
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FILE: "FILE CASE9 A 


OPTIMIZATION NO.9 
OPTIMIZATION SOLUTION 


A) PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 30000000.0 
VERT. SPAN ; 18.475 YIELD STRENGTH: 52000.0 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: 12 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
EWNEETEIY So's o's 0 0 6.6 6 640-6 6 410 4.6.6 66 0.0000 
Pe Meeetey eo c's 66. 6 oiei ee brehe ies 60.0 & 60 * 0.0000 
Sa EMPEs) OWENG co 65s) o17e 016 oo a! Oieke. 6 ose 026 0.0000 
FA... 00 Picic.s ¢/0.0.0) 0) sieice.es cio ss ~100.0000 
CONCENTRATED LOAD AT NODE... 7 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE AREA 
1 0.00000 0.00000 0.15003 0.22504 
2 2.66667 2.28090 0.15000 0.22500 
3 5.33333 4.52538 0.15025 0.22537 
& 8.00000 6.57364 0.15012 0.22518 
5 10.66666 8.44217 0.15017 0.22526 
6 13.33333 10.14746 0.15014 0.22521 
7 16.00000 11.70146 0.15023 0.22535 
8 18.66666 13.11621 0.15000 0.22500 
9 2lsisosoo 14.3974G 0.15005 0.22507 
10 23.99998 15.55445 0.15024 0.22535 
11 26.66666 16.59473 0.15000 0.22500 
l2 29.33331 17.52634 0.15000 0.22500 
is 31.99998 18.47499 0.15000 0.22500 
E)} OBJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 8.390676 
NODE NORMAL STRESS BENDING STRESS TOTAL 
Z 33701.0 0.1 33701.1 
2 33704.7 18242.0 51946.6 
3 33646.0 18393.0 52039.0 
4 33675.5 18362.1 52037.6 
5 33663.0 18380.3 52043.3 
6 33671.2 18369.3 52040.5 
7 33649.9 18380.3 52030.2 
8 33702.3 18278.8 51981.2 
9 33692.0 18333.3 52025.3 
10 33649.5 18317.5 51967.0 
11 33703.4 18112.9 51816.3 
12 33701.9 17257.1 50959.1 
13 33699.5 0.1 33699.6 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
1 1 1 0 
13 1 1 0 
G) SOLUTION VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 
1 0.000000E+00 0.000000E+00 -0.137583E-01 
2 0.273039E-01 -0.379869E-01 -0.123359E-01 
3 0.488191E-01 -0.696248E-01 -0.949834GE-02 
G 0.624682E-01 -0.935884E-01 -0.675184E-02 
5 0.696080E-01 -0.110146E+00 -0.409314GE-02 
6 0.713910E-01 -0.119528E+00 -0.150820E-02 
7 0.687846E-01 -0.121933E+00 0.101241E-02 
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0.626172E-01 
0.536328E-01 
0.424687E-01 
0.296891E-01 
0.15809GE-01 
0. 000000E+00 


~0.117539E+00 
-0. 106513E+00 
-~0.889785E-01 
-0.650603E-01 
-~0.349536E-01 

0. 000000E+00 


0.347161E-02 
0.587863E-02 
0.824613E-02 
0.105634E-01 
0.127837E-01 
0.138690E-01 
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CASE10 A 


OPTIMIZATION NO. 10 
OPTIMIZATION SOLUTION 


A) PROBLEM PARAMETERS: 


HORIZ. SPAN: 32.000 YOUNGS MODULUS: 
VERT. SPAN : 55.426 YIELD STRENGTH: 
NO OF DESIGN VAR: 24 NO OF ELEMENTS: 
B) DERIVED CONSTANTS: 
NO OF SYSTEM NODAL POINTS... 13 
NO OF DEGREES OF FREEDOM.... 39 
HORIZ. LENGTH PER ELEMENT... 2.6667 
C) STRUCTURE LOADING: 
HME suctenNGts cs ctehals.0 6 O's 6 oe 6 6 6 4 6-6 0.0000 
AVENGERS aiel sc 6. 6 Shelele © 6 ces @eee ewe 0.0000 
FM.. @eeoeeeeeeeeseeeeeeeeee ee &¢ @ @ 0.0000 
Fees e UeWeNeteielc\ e666. 6.4.0 sceievacelace « -100.0000 
CONCENTRATED LOAD AT NODE... y 
D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 
NODE X-COORD Y-COORD BASE 
i 0.00000 0.00000 0.15000 
2 2.66667 11.92522 0.15000 
3 SR EE: 19.67612 0.15000 
& 8.00000 26.02310 0.15000 
=) 10.66666 30. 88034 0.15000 
6 13. 35355 34. 78506 0.20405 
7 16.00000 39.13298 0.15000 
8 18.66666 42.25943 0.15000 
9 21.33333 45.35509 0.15000 
10 23.99998 48.17604% 0.15000 
he 26.66666 50.91873 0.15000 
12 29.33551 53.12798 0.15000 
13 31.99998 55.42598 0.15000 
E) OBJECTIVE FUNCTION: 
TOTAL STRUCTURE VOLUME: 15.148624 
NODE NORMAL STRESS BENDING STRESS 
1 37992.2 0.1 
2 38009.5 13993.7 
S 38045.6 13552.0 
& 38084.5 9528.0 
5 38055.6 1160.7 
6 27936.6 24086.0 
7 38030.6 13224.5 
8 38061.4 13822.2 
9 38059.8 13996.9 
10 38064.1 13974.1 
vO 38067.4% 4071.1 
l2 38068.1 13946.2 
13 38069.7 0. 
F) BOUNDARY CONDITIONS: 
NODE X-DISPL Y-DISPL SLOPE 
3 1 1 0 
13 if 1 0 
G) SOLUTION VECTOR: 
NODE X-DISPL Y-DISPL SLOPE 
1 0.000000E#+00 0.000000E+00 -0.101930E-01 
2 0.103071E+00 -0.389057E-01 -0.639297E-02 
3 0.129693E+00 -0.590526E-01 -0.137548E-02 
G 0.129351E+00 -0.683836E-01 -0.759869E-03 
5S 0.133686E+00 -0.787921E-01 -0.207601E-02 
6 0.135401E+00 -0.861105E-01 0.737877E-03 
7 0.122050E4+00 -0.8434G0E-01 0.261456E-02 


Ves 


30000000.0 
52000.0 
12 


AREA 
0.22500 
0.22500 
0.22500 
0.22500 
0.22500 
0.30607 
0.22500 
0.22500 
0.22500 
0.22500 
0.22500 
0.22500 
0.22500 


38069.9 
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0.111694E+00 
0.961471E-01 
0.747049E-01 
0.485993E-01 
0.262193E-01 
0.000000E+00 


-0.823635E-01 
-0.758125E-01 
-0.623201E-01 
-0.437086E-01 
-0.235824E-01 

0.000000E+00 


0.266911E-02 
0.519496E-02 
0.760786E-02 
0.844971E-02 
0.920963E-02 
0.103006E-01 
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